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ABSTRACT 


Maintenance,  a  key  element  of  Integrated  Logistics  Support,  plays  a  very 
vital  role  throughout  an  equipment/system  plaxmed  life-cycle.  Maintenance 
costs  contribute  a  msdor  portion  of  the  life  cycle  costs  of  an  equipment  or 

*  system.  Past  historical  records  have  shown  that  the  cost  associated  with 
system  maintenance  is  immense  and  usually  takes  up  a  large  portion  of  the 

% 

annual  operating  e3q)enditure.  Besides  the  costs,  sound  maintenance  efforts 
contribute  to  better  operational  availability  and  reliability  of  a  system. 
Therefore,  the  objective  is  to  attain  the  proper  balance  of  operations  between 
performance  and  effectiveness,  and  logistics  support,  which  largely  includes 
maintenance,  spares  requirements,  and  the  available  budget.  Adequate 
maintenance  is  essential  to  ensure  the  effective  and  economical  support  of  an 
equipment  or  system.  Therefore  there  is  a  need  to  design  optimal  maintenance 
policies  to  maximize  appropriate  measures  of  system  effectiveness.  These  can 
be  either  to  minimize  operational  and  maintenance  costs,  to  improve  overall 
system  reliability  or  to  maximize  operational  availability. 

In  this  thesis,  various  maintenance  scenarios  are  examined  and  the 
corresponding  optimal  maintenance  actions  are  planned  to  take  place  at 
intervals  chosen  so  as  to  maximize  an  appropriate  measures  of  effectiveness. 
Preventive  maintenance  policies  are  also  planned  so  that  the  overall  reliability 
of  the  system  is  always  kept  above  a  specified  minimum  reliability  level,  while 

•  either  keeping  the  cost  per  unit  time  to  a  minimum  or  maximizing  the 
operational  availability,  subject  to  cost  constraints. 


Ill 


Aoosssioa  Far 

RTIS  ORAfcl 

q/ 

DTIC  TAJ3 

Q 

UnaiinounoeiS 

0 

1  J.‘l  cat,  loiv...^ 

By- _ _  _ _ 

Pi) 


ii  Sail 


Dili 


TABLE  OF  CONTENTS 


I.  INTRODUCTION .  1 

A.  MAINTENANCE  .  1 

B.  MAINTENANCE  POLICIES  .  4 

C.  LIFE  DISTRIBUTIONS  .  6 

D.  FAILURE  RATE  .  7 

E.  RELIABILITY  .  8 

F.  AVAILABILITY  .  9 

G.  MAINTENANCE  TIME  DISTRIBUTION  .  10 

1.  Active  Repair  Time  .  11 

2.  Logistics  Delay  Time  .  11 

3.  Administrative  Delay  Time  .  12 

4.  Distribution  .  12 

II.  OPTIMAL  MAINTENANCE  POLICIES  AND  MATHEMATICAL 

MODELS  .  18 

III.  POLICY  KAGE  REPLACEMENT)  .  17 

A.  NUMERICAL  ILLUSTRATION  .  23 

B.  SENSITIVITY  ANALYSIS  .  30 

1.  Alpha  .  30 

2.  Lambda . . .  33 

3.  Costs .  36 

4.  Dovmtime  . 37 

C.  EFFECTS  OF  SIMPLE  PREVENTIVE  MAINTENANCE  ...  38 

1.  Case  1  .  41 

2.  Case  2 .  43 

3.  Numerical  Illustration .  46 


4 


iv 


IV.  POLICY  II  (MINIMAL  REPAIR  WITH  AGE  REPLACEMENT)  .  .  60 


A.  COST  MODEL  .  61 

B.  AVAILABILITY  MODEL  .  69 

V.  POLICY  III  (MINIMAL  REPAIR/PAILURE  REPLACEMENT/ 

PREVENTIVE  REPLACEMENT) .  74 

A.  COST  MODEL  .  77 

1.  Case  1  .  79 

,  2.  Case  2  .  79 

3.  Weibull  Example  .  81 

«  4.  Numerical  IHustration .  Sr 

B.  AVAILABILITY  MODEL  .  86 

C.  MAXIMIZE  AVAILABILITY  SUBJECT  TO  BUDGET 

CONSTRAINT .  101 

VI.  CONCLUSIONS  AND  RECOMMENDATIONS  . 104 

APPENDKA  COMPITTER  PROGRAM  AND  OUTPUT  . 108 

APPENDIX  B  DERIVATION  OP  COST  RATE  AND 

AVAILABILITY  FUNCTIONS . 128 

APPENDIX  C  DESCRIPTION  OF  IM8L  SUBROUTINE 

DQDAG  . 133 

APPENDIX  D  DESCRIPTION  OP  IMSL  SUBROUTINE 

DCSINT  . 136 

APPENDIX  E  DESCRIPTION  OP  IMSL  FUNCTION 

DCSITG . 137 

LIST  OP  REFERENCES  . 138 

.  BIBILIOGRAPHY  . 140 

INITIAL  DISTRIBUTION  LIST  . 141 

t 


V 


ACKNOWLEGEMENT 


I  am  most  grateAil  to  my  thesiB  advisor*  Professor  Donald  P.  Qaver,  for 
his  sound  and  professional  advice,  guidance,  assistance  and  encouragement  in 
the  completion  of  this  thesis. 

I  would  also  like  to  express  my  appreciation  to  my  second  reader, 
Professor  Lyn  Whitaker  for  her  good  comments,  and  also  to  Professor  Patricia 
Jacobs  for  her  assistance. 

Lastly  I  am  very  grateful  to  my  lovely  wife,  Telly,  for  her  support  and 
patience,  and  also  to  my  two  lovely  sons  Manharpreet  and  Sandeep. 


vi 


THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed  in  this 
research  may  not  have  been  exercised  for  all  cases  of  interest.  While  every 
effort  has  been  made,  within  the  time  available,  to  ensure  that  the  programs 
are  free  of  computational  and  logic  errors,  they  cannot  be  considered  validated. 
Any  application  of  these  programs  without  additional  verification  is  at  the  risk 
of  the  user. 
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L  INTRODUCTION 


A.  MAINTENANCE 

Maintenance  constitutes  a  series  of  actions  to  be  taken  to  restore  or  retain 
an  equipment  or  system  in  an  effective  operational  state.  [Hef.  1]. 

Maintenance  occurs  at  three  levels: 

•  Oiiganizational  Maintenance. 

•  Intermediate  Maintenance. 

•  Depot  Maintenance. 

The  Organizational  Maintenance  activity  is  performed  by  the  ship’s  staff, 
for  frigates  or  equivalent,  using  the  onboard  spares,  tools,  test  equipment,  and 
the  documentation  furnished.  For  smaller  ships  and  patrol  craft. 
Organizational  Maintenance  is  carried  out  with  the  assistance  of  the  Squadron 
Support  Teams.  Intermediate  Maintenance  is  performed  at  fleet  workshops, 
while  Depot  Maintenance  is  carried  out  at  the  naval  dockyard  or  by 
commercial  shipyards  and  firms. 

Any  equipment  or  system  introduced  into  the  service  needs  to  be 
maintained  adequately  so  that  it  is  readily  available  to  perform  a  mission 
successfully  at  an  acceptable  performance  level.  Besides  that,  maintenance 
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helps  to  extend  the  useful  life  of  the  eqmpment/system  and  ensures  the  safety 
of  personnel  using  it. 

Basically,  maintenance  can  be  broadly  divided  into  two  types: 

*  Preventive  Maintenance. 

*  Corrective  Maintenance. 

Figure  1.1  shows  the  various  maintenance  forms. 


Figure  1.1  The  Forma  of  Maintenance 


Preventive  maintenance  is  scheduled  maintenance  that  is  organized  and 
carried  out  in  accordance  with  the  documentation  and  records  at  a 
predetermined  time  following  a  predetermined  plan.  Preventive  maintenance 
is  normally  associated  with  maintenance  that  oc(mrs  when  the  system  is  still 
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operating.  The  main  aim  of  preventive  maintenance  actions  is  to  improve  the 
overall  system  reliability  and  to  avoid  sudden  or  unexpected  failures  which 
may  be  cata  strophic.  Preventive  maintenance  can  be  sub-divided  into: 


■  Simple  Preventive  Maintenance. 

■  Complete  Overhaul/Preventive  Replacement. 

Simple  preventive  maintenance  actions  usually  consist  of  inspections, 
adjustments,  tuning,  cleaning,  lubrication,  minor  calibration,  and  replacement 
of  worn  out  components  and  parts  before  they  actually  fail.  This  type  of 
maintenance  is  usually  performed  at  the  Organizational  Maintenance  level  and 
often  it  does  not  affect  the  downtime  of  the  equipment/system.CThe  downtime 
for  these  maintenance  actions  is  usually  negligible). 

Complete  overhauls  are  carried  out  to  bring  the  state  of  the  deteriorating 
system  back  to  "as  good  as  new  condition",  and  to  prevent  mqjor  impending 
failures.  This  type  of  maintenance  is  usually  carried  out  at  the  Depot 
Maintenance  level  during  planned  refit  periods  and  it  usually  incurs  some 
significant  downtime.  Preventive  replacements  are  carried  out  when  the 
overhaul  is  not  economical  or  when  the  parts  are  non-repairable. 

Con'ective  maintenance  is  the  unscheduled  maintenance  carried  out  to 
restore  a  failed  system  to  its  operating  state.  Corrective  maintenance  can  be 
subdivided  into; 
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*  Minimal  Repair. 

•  Mc^or  Overhaul/Failure  Replacement. 

Minimal  repairs  are  minor  repairs  or  component  replacements  carried  out 
on  failed  components  or  assemblies  that  restore  the  system  to  its  operational 
state  without  significantly  improving  the  overall  condition  of  the  system.  This 
type  of  maintenance  is  usually  carried  out  at  the  Organizational  Maintenance 
level, 

Mc^or  overhauls  are  carried  out  when  the  system  experiences  a  sudden 
mayor  failiure  and  the  work  has  to  be  carried  out  at  the  Intermediate 
Maintenance  level  or  at  the  Depot  Maintenance  level.  Failure  replacement 
usually  occurs  when  the  system  is  beyond  economical  repair  or  when  the  parts 
are  non*repairable. 

B.  MAINTENANCE  POLICIES 

Our  main  goal  is  to  design  a  maintenance  policy  that  will  maximize 
appropriate  measures  of  system  effectiveness;  this  can  be  either  to  minimize 
operational  and  maintenance  costs,  to  improve  overall  system  reliability  or  to 
maximize  operational  availability  (uptime).  The  most  relevant  measure  of 
effectiveness  used  is  to  maximize  a  suitable  measure  of  operational  availability 
subject  to  budget  constraints.  Ideally,  we  would  like  to  carry  out  many 
preventive  maintenance  and  inspection  routines  to  ensure  that  the  system  is 
in  an  optimum  operating  condition,  to  avoid  sudden  and  catastrophic  failures 
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and  damage  to  the  system,  and  also  to  prevent  accidents  which  may  be 
detrimental  to  the  people  working  in  the  vicinity.  On  the  other  hand,  excessive 
preventive  maintenance  actions  may  not  only  be  unnecessarily  costly,  requiring 
many  manhours  but  may  actually  prematurely  age  equipment.  Therefore,  we 
need  a  balance  or  trade-off  between  two  extremes.  Our  objective,  then,  is  to 
select  an  optimal  maintenance  policy  for  the  particular  equipment  or  system, 
and  to  decide  when  to  carry  out  the  associated  maintenance  routines.  Some 
conceivable  maintenance  policies  are  the  following: 


*  No  preventive  maintenance  and  corrective  maintenance  is  carried  out. 
The  system  is  replaced  at  a  fixed  age.  When  the  system  fails  before  this 
age  it  is  replaced  with  a  new  one.  This  type  of  maintenance  is  normally 
associated  with  low-level  subsystems:  items  such  as  components,  sealed 
modules  or  other  non-repairable  parts.  Some  common  examples  are  the 
magnetron  in  the  radar  transmitter  unit  and  the  belts  found  in  motors. 

*  No  preventive  maintenance  is  carried  out.  The  system  is  replaced  at  a 
fixed  age.  For  failures  that  occur  in  between  the  planned  replacement 
age,  minimal  repairs  are  carried  out  to  restore  the  system  to  an 
operational  state.  This  type  of  maintenance  is  usually  associated  with  a 
system  consisting  of  several  components  such  that  when  a  component 
fails  the  system  fails,  and  replacing  the  component  restores  the  system 
to  operation.  Minimal  repairs  are  often  carried  out  in  the  field,  i.e. 
onboard  ships  or  at  a  forward  air  base. 

*  The  system  is  renewed  at  a  fixed  age  either  by  replacement  or  overhaul 
after  which  the  system  is  "as  i^od  as  new".  Failures  in  between  the 
planned  age  replacement  can  be  classified  into  type  I  and  type  II.  Type 
I  failures  are  simple  failures  which  are  remedied  by  minimal  repairs 
using  the  support  elements  onboard  ships,  and  II  failures  are  mtgor 
failures  which  require  base  support  faciUties  and  are  rectified  by  part 
replacement  or  complete  overhaul. 

*  Both  preventive  and  corrective  maintenance  actions  are  employed. 
Preventive  maintenance  is  planned  at  some  time  interval  to  improve  the 
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,  Other  classes  of  distributions  which  exhibit  wear  out  or  wear-in  such  as 
NBU  (new  bettter  than  used)  or  UBN  (used  bettter  than  new)  can  also  be  used 
and  are  discussed  in  detail  in  Barlow  and  Proscham.  [Ref.  2]. 

Most  of  the  systems  which  are  newly  installed  onboard  the  ships  initially 
have  a  decreasing  failure  rate,  sometimes  referred  to  as  the  infant  mortality 
phase  or  running-in  period.  Then  the  failure  rate  becomes  rather  constant  for 
some  time  and  finally  increasing,  exhibiting  wear  out.  In  reliability,  such 
failure  rate  functions  are  said  to  have  a  "bathtub"  shape.  The  running-in 
period  is  usually  under  contractual  obligation. 

D.  FAILURE  RATE 

The  failure  rate  or  hazard  rate  is  one  of  the  most  important  statistical 
characteristics  of  any  equipmenVsystem  frequently  used  in  maintenance  or 
replacement  studies. 

The  failure  rate,  h(t),  is  usually  defined  as  the  "instantaneous"  conditional 
probability  of  failure  at  age  t,  given  that  it  has  survived  to  age  t.  When  F  has 
density  f,  h(t)  is  given  by 

f(t) 

h<t)  -  -  (1.1) 

1  -  F(t) 

The  failure  rate  h(t)  can  also  be  expressed  in  terms  of  the  hazard 
function,  H(t) 
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dH(t) 

h<t)  - -  (1.2) 

dt 

The  hazard  function  is  related  to  the  distribution  function  by  the 
following  relationship 

H<t)  -  -  log  {I  -  F(t)}  (1.8) 


E.  RELIABILITY 

When  deciding  upon  system  maintenance  polides,  the  frequency  of 
maintenance  actions  becomes  a  significant  parameter.  The  frequency  of 
maintenance  for  a  given  system  is  highly  dependent  on  the  reliability  of  that 
system.  In  general,  as  the  reliability  of  a  system  increases,  the  frequency  of 
maintenance  actions  will  decrease;  conversely,  the  frequency  of  maintenance 
actions  will  increase  as  system  reliability  is  degraded.  [Ref.  1] 

The  reliability  function  or  survival  function  R(t)  is  given  by 

JKO  ■  I  -  W -/'AO*  (M) 

R(t)  is  the  probability  that  a  new  system  will  perform  its  mission 
satisfactorily  for  at  least  a  certain  time  t.  If  T  is  the  time  to  failure  of  the 
system  then 

R(t)  -  P(T  >  t)  (1.5) 

Sometimes  we  are  interested  in  the  chance  of  survival  of  the  system  in 
the  future  given  that  it  has  survived  up  till  now.  This  is  called  the  conditional 
survival  function.  So  if  T  is  the  time  to  failtire  of  the  system  with  the  survival 
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function  [1  -  F(t)],  then  the  reliability  of  the  system,  R(t;a),  at  some  age  a  S  0 
is  given  by 


Cl  -  r(t+a)l 

■  -  (1.6) 

[1  -  r(«)l 


F.  AVAILABILITY 

Availability  is  a  measure  of  system  readiness  and  it  is  one  of  the  most 
important  measures  of  effectiveness  usually  employed  in  mission-oriented 
situations  especially  in  the  military  environment.  Operational  availabiUty  is 
the  probability  that  a  system  or  equipment,  when  used  under  stated  conditions 
in  an  actual  operational  environment,  will  operate  satisfactorily  when  called 
upon  [Ref  1].  System  availability  is  influenced  both  by  the  inherent  failure- 
proneness  of  the  system  and  by  the  time  and  resources  (support  elements)  it 
takes  to  restore  a  failed  system  to  service.  [Ref.  3].  Times  to  failure  or  ’up 
times’  and  to  restoration  or  ’down  times’  may  vary  considerably,  and  not 
necessarily  independently,  depending  upon  the  mode  of  failure,  the  time 
required  to  diagnose  the  failure,  availability  of  special  tools,  test  equipment, 
and  spare  parts,  and  the  proper  dcctucnentation  and  the  required  personnel 
skills.  The  long-run  availability  or  steady  state  is  expressed  as  follows: 

1[U] 

A,  - -  (1.7) 

E[01  +  l[D3 
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where 


ECU]  is  the  expected  uptime  of  the  system 
E[D]  is  the  expected  downtime  of  the  system 

For  a  system  operating  at  sea  the  availability  at  equation  (1.6)  is  usually 
expressed  as 

MTBM 

K  - - -  (1.8) 

MTBM  .f  MDT 

where 

MTBM  is  the  mean  time  between  maintenance 
MDT  is  the  mean  downtime,  which  includes  the  mean  active 
maintenance  time  (M),  expected  logistics  delay  time  (LDT)  and  the  expected 
administrative  delay  time  (ADT).  The  mean  active  maintenance  time  includes 
the  expected  time  for  preventive  and  corrective  maintenance. 

G.  MAINTENANCE  TIME  DISTRIBUTION 

The  time  required  to  carry  out  simple  preventive  maintenance  actions  and 
overhauls  can  generally  be  modeled  as  normally  distributed  with  mean  p  and 
standard  deviation  o,  Most  of  tliese  tasks  are  standard  and  are  carried  out  in 
accordance  with  the  planned  maintenance  schedules;  which  stipulates  the 
procedures  to  follow,  spares,  material,tools  and  test  equipment  that  are 
required  to  perform  the  maintenance  actions.  [Ref  4].  The  tasks  usually 
requires  a  fixed  sunount  of  time  to  accomplish  with  very  little  variation. 
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The  time  required  for  corrective  maintenance  actions  can  be  divided  into 
three  basic  categories: 

*  Active  repair  time. 

*  Logistics  delay. 

*  Administrative  delay. 

1.  Active  Repair  Time 

Active  repair  time  depends  on  the  environment^  state  of  equipment 
(hot  or  cold),  and  skill  level  of  the  technician;  and  it  can  be  sub-divided  into  the 
following  categories: 

*  Recognition  or  detection  (  often  the  time  until  actual  occurrence  of  a 
failure  and  its  recognition  is  not  known). 

*  Localization  or  diagnosis. 

*  Correction  or  repair. 

«  Verification  or  check. 

2.  Logistios  Delay  Time 

Logistics  delay  time  constitutes  downtime  that  is  expended  while 
waiting  for  the  availability  of  a  spare  part,  waiting  for  a  special  tool  or  test 
equipment  to  perform  repair,  waiting  for  transportation,  and  waiting  to  use  a 
facility  required  for  the  repair. 
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a.  Administrative  Delay  Time 

Administrative  delay  time  constitutes  downtime  of  adminstrative 
nature,  such  as  personnel  assignment  priority  and  organizational  constraint. 

4.  Distribution 

The  distributions  most  commonly  used  to  describe  the  downtime  for 
corrective  maintenance  actions  are  exponential  and  log-normal.  The 
exponential  distribution  tends  to  fit  the  type  of  equipment  that  requires 
relatively  short  durations  of  repair  and  usually  corresponds  to  the  replacement 
of  a  failed  unit.  Occasionally,  much  longer  times  may  be  required  for  nu^or 
repair  or  for  spares.  The  lognormal  distribution  is  usefiil  for  situations  where 
there  are  few  downtimes  of  short  duration,  a  large  number  of  downtimes 
closely  grouped  about  some  modal  value  and  a  few  downtimes  of  long 
durations. 

If  X,  the  downtime,  is  a  random  variable  having  the  lognormal 
distribution  given  by  the  probability  density  ftinction 

..flajt  - 

/*(x)  ■  — =r*  0  s  X  <  (1.9) 

x/2no’ 

then  the  logarithm  of  X  is  Normal  with  mean  p  and  variance  cr*.  In  the  thesis 
We  shall  study  several  maintenance  policies  applicable  to  systems  onboard 
ships.  In  Chapter  II  we  introduce  optimal  maintenance  for  three  types  of 
policies.  These  are  considered  in  detail  in  Chapters  III,  IV  and  V. 
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n.  OPTIMAL  MAINTENANCE  POUCIES  AND  MATHEMATICAL 

MODELS 

In  the  commercial  environment  determination  of  the  optimum 
maintenance  policy  and  time  is  of  great  economic  importance.  However,  in 
many  military  sitiiations,  failure  of  a  system  in  an  operational  environment  is 
not  only  going  to  be  more  costly,  but  dangerous  and  may  jeopardise  the  success 
of  a  mission.  If  a  system  has  an  increasing  Mure  rate,  such  as  the  failure  rate 
of  a  WeibuU  distribution  with  shape  parameter  >  1,  it  may  be  wise  to  replace 
or  overhaul  the  system  before  it  has  aged  too  greatly.  [Ref.  6].  This  is  very  true 
for  systems  onboard  the  ships  especially  when  they  are  operating  many 
hundreds  of  miles  from  their  home  base  and  hence  the  support  elements  are 
not  completely  and  readily  available.  A  Mure  at  sea  may  be  catastrophic,  in 
terms  of  cost  and  operational  requirements.  Although  we  can  not  completely 
avoid  fkilures,  however  we  cem  reduce  the  chance  or  probability  of  such 
catastrophic  failures.  This  can  be  done  by  studying  the  Mure  distributions 
and  then  employing  appropriate  maintenance  actions  to  maximize  the  various 
measures  of  effectiveness. 

The  appropriate  maintenance  actions  could  be  either  preventive 
maintenance,  failure  replacement,  minimal  repair  and  preventive  replacement 
or  complete  overhaul.  The  objective  of  a  maintenance  policy  is  to  find  a 
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sequence  of  times  for  carr^dng  out  the  various  maintenance  actions  that 
maximizes  the  appropriate  measures  of  effectiveneBS  over  the  operational  and 
maintenance  cycle  of  the  system. 

The  most  commonly  used  maintenance  policy  is  the  policy  based  on  age, 
usually  referred  as  the  age  replacement  polity.  Sometimes  the  maintenance 
policy  is  based  on  the  running  hours  of  the  system  or  equipment.  If  a  system 
consists  of  many  identical  components,  then  the  maintenance  of  these 
components  are  done  in  a  block  or  group  and  is  called  the  block  replacement 
policy,  such  as  the  replacement  of  the  diodes  in  the  exciter  unit  of  an 
alternator  where  the  accessibility  is  poor. 

A  recent  survey  of  Preventive  Maintenance  Models  for  stochaiticBlly 
deteriorating  single>unit  system  [Ref.  6]  highlighted  the  use  of  some 
optimization  models  for  repair  and  replacement  policy  evaluation.  Most  of 
these  models  were  based  on  minimizing  the  long-run  expected  costs  per  unit 
time  of  replacement  and  minimal  repair  as  the  measure  of  effectiveness.  The 
basic  minimal  repair  model  developed  by  Barlow  and  Hunter  [Ref.  7]  has  been 
generalised  and  modified  by  many  authors  to  fit  more  realistic  situations. 

Minimal  repair  models  generally  assume  the  following;  [Ref.  6]. 

•  The  system’s  failure  rate  fimction  is  increasing. 

•  Minimal  repairs  do  not  affect  the  failure  rate  of  the  system. 
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•  The  cost  of  a  minimal  repair  Cf  is  less  than  the  cost  of  replacing  the 
entire  system  Cp. 

•  System  failures  are  immediately  detected. 

The  long*run  expected  cost  per  unit  time  using  a  replacement  age  t  for 
the  basic  model  is  given  by 

.  (2.1) 

t 

where  N(t)  represents  the  expected  number  of  failures  (minimal  repairs) 
during  the  period  (0,t]. 

Using  the  basic  minimal  repair  model  as  developed  by  Barlow  and  Hunter 
[Ref.  7],  Tilquin  and  Cldroux  [Ref.  8]  investigated  an  optimal  replacement 
policy  for  the  case  where  an  adjustment  cost  C.(ik),  incurred  at  age  ik, 
i  ■  1 ,2,3,...  and  k  >  0,  is  added  to  the  basic  costs  Cp  and  Cf.  They  showed  that 
the  long*run  expected  cost  per  unit  time  is  given  by 

.  c,m  *c,  *  c-j.m 

t 


where 


and  v(t)  rfipresents  the  number  of  adjustments  in  the  period  (0,t].  Tilquin  and 
Cl^roux  [Ref.  8]  showed  that  the  global  minimum  for  equation  (2.2)  exists  in 
the  interval  [0,m]  when  the  life  distribution  is  IFR. 
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In  the  thesis  we  will  examine  the  various  maintenance  scenarios  here 
culled  the  policies  and  then  based  on  these  we  will  formulate  the  appropriate 
mathematical  models  using  the  stochastic  and  reliability  theory.  This 
mathematical  modek  will  depend  on  the  desired  measures  of  efifectiveness 
required.  From  these  models  we  can  obtain  the  times  for  carrying  out  the 
appropriate  maintenance  actions.  The  measure  of  effectiveness  that  will  be 
considered  are: 

*  Minimizing  the  costs. 

*  Maximizing  availability. 

*  Mission  reliability. 

The  maintenance  actions  will  be  different  for  various  measures  of 
effectiveness  and  it  is  up  to  the  Decision  Maker  to  select  which  one  is  suitable 
for  his  scenario. 

The  following  maintenance  policies  are  of  interest; 

*  Policy  I  Age  Replacement. 

*  Policy  II  Minimal  Repair  with  Age  Replacement. 

*  Policy  III  Minimal  Repair,  Failure  Replacement  / 

Overhaul,  Preventive  Replacement  / 

Overhaul. 

These  policies  are  discussed  in  the  following  chapters. 
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HL  POUCY I  (AGE  REPLACEMENT) 


The  system  is  replaced  at  the  time  of  failure  or  at  some  fixed  time  tp 
whichever  comes  first.  (Instead  of  replacement,  we  could  also  overhaul  the 
system, which  on  completion  of  overhaul  is  assumed  as  "good  as  new"). 

This  type  of  scenario  is  usually  associated  with  the  repair  by  replacement 
policy  often  adopted  at  sea  and  applies  to  modules  and  sub-assemblies  of 
equipment  which  requires  support  elements  not  available  at  sea.  In  lieu  of 
repairing  these  modules  and  sub-assemblies  at  sea,  sufficient  spares  are 
carried  onboard  or  propositioned  at  the  forward  operating  areas  so  as  to 
accomplish  the  respective  missions  suooessfhlly.  The  optimal  number  of  spare 
requirements  are  based  on  the  measures  of  effectiveness  desired. 

The  age  at  which  the  operating  system  is  replaced  depends  on  the 
following  factors: 

*  Failure  distribution. 

*  Costs  of  failure  and  preventive  replacement. 

*  Downtime  of  failure  and  preventive  replacement. 

*  Measure  of  effectiveness: 

•  Minimise  costs. 

•  Maximize  availability. 

•  Mission  reliability. 
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When  evaluating  this  policy  the  following  assumptions  are  made; 


•  Planned  replacements  are  less  costly  than  failure  replacements. 

•  The  mean  downtime  for  a  planned  replacement  is  less  than  that  for  a 
fa^ure  replacement.  An  unexpected  failure  may  incur  additional  Liogistics 
Delay  Time  (LDT)  and  Administrative  Delay  Time  (ADT),  especially  if  the 
fedlures  occured  whilst  the  ship  is  at  sea. 

•  The  system  exhibits  an  increasing  failure  rate  distribution  i.e.  h(t) 
increases  as  t  increases. 

•  The  cost  and  downtime  associated  with  simple  preventive  maintenance 
actions  and  minimal  repairs  is  negligible, 

•  Preventive  maintenance  actions  and  minimal  repairs  do  not  improve  the 
reliability  of  the  system. 


Let  T  be  the  time  to  fkilure  of  the  system.  T  is  a  random  variable  with  an 
IFR  distribution  Amotion  F(t).  A  cycle  is  completed  everytime  a  replacement 
is  made.  It  can  be  either  a  failure  replacement  or  a  preventive  replacement. 
(The  system  probabilistically  starts  over  again  and  each  replacement 
constitutes  a  renewal).  By  using  the  Renewal  Reward  Process,  Ross  [Ref.  9]: 
Expected  long  run  average  cost,  C(tp),  is  given  by 

Kxpaotad  cost  laourrad  during  a  oyela 

C  (tp)  «  - - - 

Ixpaotad  length  of  oyole 


l[L] 
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where 


C  is  the  total  maintenance  cost  incurred  during  a  cycle 
L  is  the  lengtli  of  a  cycle 
Let  tp  be  the  planned  replacement  age 


failure 

replaoeaent 


preventive 

replacement 


—T—"— 


-tp- - 


failure 

replacement 


C 


(  Cp 
{ 

{  c, 


T  >  tp 

T  S  tp 


(8.2) 


where  Cp  is  the  cost  of  preventive  replacement 
Cf  is  the  cost  of  failure  replacement 
Cr  >  Cp 

{  Rp  +  tp  T  >  tp 

L  -  { 

(  R,  +  T  T  i  tp 

where  Rp  is  the  time  of  preventive  replacement 
Rf  is  the  time  of  failure  replacement 
Rf  >  Rp 

Then, 


C(fp) 


Cp[i-i»i(>p)l^qF(tp)] 


(8.8) 


(3.4) 
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and  the  above  equation  can  be  simplified  to  the  standard  form  as  shown  in 
Barlow  [Ref.  7]  and  Jardine  [Ref.  10].  See  Appendix  B  for  details. 


where  M(tp)  is  the  mean  life  during  a  cycle  and  is  given  by 


(3.5) 


«(V  -  0js> 

The  optimal  preventive  replacement  age  can  be  fmmd  by  finding  the  value 
of  tp  that  minimizes  the  cost  function  in  equation  (3.4)  and  the  optimal  value 
of  t*p  is  that  value  that  satisfies  the  following  equation  [Appendix  B] 

^  *  *(ip  m 

The  other  measure  of  effectiveness  that  the  Decision  Maker  is  often 
interested  is  the  aiailability  (the  probability  the  system  is  up  at  any  time  t). 

Assuming  that  there  are  only  two  states,  that  is  the  system  is  either  up 

or  down,  then  from  the  Regenerative  Process  Ross  [Ref  9],  Availability ,A(t),  is 

the  Expected  amount  of  time  the  system  is  up  during  a  cycle  divided  by  the 

Expected  time  of  a  cycle  and  is  given  by 

Neiua  life  during  a  eyela 

A(t)  - - 

Knpeeted  length  of  a  cycle 
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m;> 


(3.8) 


A(t;> 


The  optimal  preventive  replacement  age  can  be  found  by  finding  the  value 
of  tp  that  maximizes  the  availability  function  in  equation  (3.8)  and  the  optimal 
value  of  t*p  is  that  value  that  satisfies  the  following  equation.  [Appendix  B] 


m;>  -  [ 


(3.9) 


The  other  measure  of  effectiveness  of  interest  is  the  mission  reliability 
which  is  defined  to  be  the  probability  that  the  system  will  complete  a  certain 


mission  of  duration  d  when  it  is  at  age  t.  This  is  given  as  follows 

R(t  4  d) 

R(t,d)  -  - 

R(t) 

1  -  r(t  4  d) 

■  -  (8,10) 

1  -  r(fc) 


In  military  applications  we  usually  like  to  maximize  the  availability  of  the 
system  subject  to  some  budget  constraint  C: 
siax  A(tp) 

a.t.  C(tp)  ^  C  (8.11) 

tp  >  0 

Some  of  the  continuous  life  distributions  that  are  commonly  used  to  model 
the  increasing  failure  rate  of  the  system  are: 
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a.  Weibull. 


b.  Gamma  with  shape  parameter  >  1. 

c.  Log-Normal  (depends  on  parameter). 

If  we  do  not  know  the  probability  life  distribution  then  we  have  to  resort 
to  a  nonparametric  approach  [Ref.  10]. 

It  should  be  noted  that  if  the  failure  rate  of  the  system  is  constant  i.e.  it 
exhibits  ihe  Exponential  distribution  then  by  virtue  of  its  memoryless  property 
we  do  not  cany  out  any  preventive  replacement  no  matter  what  because  a  new 
system  is  just  as  bad  and  good  as  an  old  one. 

For  our  ease  let  us  assume  that  from  historical  data  we  know  that  this 
system  has  a  Weibull  distribution  with  shape  parameter  u  and  scale  parameter 
X.  Then  the  probability  distribution  function  F(t)  is  given  by 

a 

-(Xr) 

FCO-W  a>l  »X>0  .TkO  (3.12) 

and  the  failure  rate  h(t)  is  given  by 

/!i(0"aX  t  a>l  ,X>0  ,rkO  (3.13) 

and  the  hazard  function  H(t)  is  given  by 

a 

/r(0  -  (xr)  a>u  x>0,  uo  (3.14) 
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It  should  be  noted  that  if  the  shape  parameter  alpha  »  1.0,  then  the 
failure  rate,  h(t)  b  "K,  which  is  a  constant  and  this  corresponds  to  an 
Exponential  distribution  whose  distribution  hmction  is  given  by 

-  1  -  X>0/a0  (3.15) 

Figures  3.1  and  3.2  shows  the  plot  of  the  Weibull  density  function  and  the 

corresponding  hazard  function  respectively.  It  can  be  observed  that  as  alpha 
increases  the  failure  rate  increases  more  rapidly. 

It  should  be  noted  that  for  a  Weibull  distribution 

a 

has  no  simple  elementary  closed  form  solution,  but  can  be  expressed  in  terms 
of  the  incomplete  Gamma  fhnction. 

A.  NmiERICAL  ILLUSTRATION 

A  case  example  was  taken  with  the  parameters  as  follow: 

Cp  B  $  26000 

CrB$  37600 
Rp  B  8  hours 
Rf  B  16  hours 
a  ■  3.0 

X  B  1/1390  hours 
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Figure  8.2  Hazard  Function  of  the  Welbull  dlitribution  H(t),  a  >  1 

and  1/lambda  «  1890  bourn 
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The  above  example  was  extracted  from  [Ref.  12].  This  example  was  based 
on  the  overhaul  of  aircraft  engines.  This  same  principle  can  be  applied  to  the 
overhaul  of  motors,  compressors  pumps  and  weapon  systems  onboard  ships. 

An  1M8L  subroutine  DQDAO  was  available  on  the  1BM3033P  main  frame 
computer  at  the  Naval  Postgraduate  School,  [Ref.  13],  which  could  integrate 
the  fimction  very  accurately  and  expeditiously  by  using  a  globally  adaptive 
scheme  based  on  Qauaa-Kronrod  rules.  The  estimate  of  the  absolute  value  of 
the  error  using  this  scheme  was  10*^^  A  brief  description  of  this  subroutine  is 
shown  in  Appendix  C. 

The  cost  ftinction  and  the  availability  function  as  in  equations  (3.6)  and 
(3.8)  respectively  were  plotted  against  time.  Figures  3.3  and  3.4  shows  these 
plots  respectively.  It  is  observed  that  the  cost  is  minimum  at  about  tp*  ■  1463 
hours  and  at  a  rate  of  about  $28.96  per  hour.  At  this  critical  value  of  tp  (1463 
hours)  the  availability  is  about  0.9883.  It  is  also  observed  that  the  availability 
is  maximum  at  tp'  ■  1126  hours  giving  an  availability  of  0.9888  .  For  this 
maximum  value  of  availability  we  need  a  budget  of  at  least  $29.92  per  hour. 
Based  on  this  information  it  is  up  to  the  Decision  Maker  to  choose  the  optimal 
time  tp*  to  carry  out  preventive  replacement  or  overhaul. 

From  the  graphs  in  Figures  3.3  and  3.4  it  is  observed  that  the  plots  are 
fairly  flat  and  this  gives  some  flexibility  to  the  Decision  Maker.  The 
replacement  age  need  not  occur  exactly  at  the  optimal  age. 
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Figure  8.8  Plot  of  the  Cost  Rate  Funotton*  C(tp). 


Figure  8.4  Plot  of  the  Availability  Function,  ACt^). 
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The  optimal  values  of  tp*  for  the  minimum  cost  and  the  maximum 
availability  as  mentioned  above  also  satifies  the  equations  (3.7)  and  (3.9) 
respectively.  The  values  are  1463.46  and  1126.38  hours  respectively. 

The  Decision  Maker  may  also  be  interested  in  the  mission  reliability  of 
the  system  /  equipment.  The  graph  in  Figure  3.6  shows  the  probability  that  the 
system  /  equipment  will  sustain  24  hours  of  continuous  operation  successfully 
when  the  equipment  is  at  a  certain  age,  §ay  t  hours.  Suppose  the  Decision 
Maker  wants  a  reliability  of  not  less  than  0.96;  then  the  optimum  replacement 
age  t^*  is  1371  hours. 


Figure  8.6  Plot  of  the  Reliability  Function,  R(t). 


From  the  plot  in  Figure  3.6  it  is  observed  that  the  reliability  decreases  as 
the  equipment  ages.  This  shows  that  even  if  the  equipment  does  not  fail  and 
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if  it  has  to  sustain  an  operation  for  a  specified  duration  of  time,  an  older 
equipment  ^11  be  less  reliable  and  thus  may  not  fbUlll  an  operational  mission 
successfully. 

Table  3.1  shows  the  values  of  the  cost  rate,  availability  and  the  reliability 
for  the  replacement  age  around  tho  optimal  values.  It  is  observed  that  the 
availability  value  is  very  close  whether  the  replacement  age  is  at  900  hours  or 
at  1600  hours.  This  is  mainly  attributable  to  the  small  values  of  the  downtimes 
for  the  preventive  replacement,  Rp,  and  Mure  replacement,  R,. 

The  three  measures  of  effectiveness  as  discussed  above  are  important 
criteria  to  determine  the  amount  of  spares  required  to  be  carried  onboard  or 
prepositioned  at  forward  operating  areas.  When  a  ship  is  assigned  to  an 
operational  area,  say  for  2  months,  away  from  the  homeport,  the  ship  will  have 
to  rely  on  the  forward  base  for  replenishment.  The  replenishment  cycle  occurs 
every  two  or  three  weeks  and  usually  takes  about  two  days.  So  if  the  optimal 
replacement  age  is  close  (not  necessarily  within)  to  these  stand-off  periods  we 
can  undertake  the  replacement  actions  during  these  periods.  It  should  be  noted 
that  if  we  do  not  carry  the  replacement  action  when  it  is  due  then  if  we 
continue  to  delay  these  actions  everytime  they  are  due,  on  the  long  run  we  are 
going  to  experience  high  costs,  poor  availability  and  low  reliability. 

Since  we  have  three  measures  of  effectiveness  and  the  respective  optimal 
replacement  ages,  it  is  up  to  the  Decision  Maker  to  decide  which  one  he  is 
going  to  give  top  priority. 
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TABLE  3.1  VALUES  OF  COST/HR,  AVAILABILITY  AND 
RELIABlLlTy  FOR  VALUES  OF  REPLACEMENT  AGE  CLOSE  TO 

THE  OPTIMUM. 


Replacement 


Coet 


Age  (tp*) 
900.0 
960.0 
1000.0 
1060.0 
1100.0 
1160.0 
1200.0 
1260.0 
1300.0 
1360.0 
1400.0 
1460.0 
1600.0 
1660.0 
1600.0 


($/hr) 

32.78 

31.90 

31.18 

30.69 

30.12 
29.76 
29.47 
29.26 

29.12 
29.02 
28.97 
28.96 
28.96 
28.99 
29.04 


Availability 

(A) 

0.9884 


0.9886 

0.9887 

0.9888 

0.9888 

0.9888 

0.9888 

0.9887 

0.9886 


0.9886 


0.6886 

0.9884 

0.9883 

0.9882 

0.9881 


ReliabiUty 

(R) 

0.9779 

0.9766 

0.9729 

0.9702 

0.9674 

0.9644 

0.9614 

0.9682 

0.9649 

0.9616 

0.9480 

0.9443 

0.9406 

0.9367 

0.9327 


29 


B.  SENSITIVITY  ANALYSIS 

Since  we  do  not  actually  know  the  fiedlure  parameters,  it  is  also  important 
to  do  sensitivity  analysis  and  see  how  the  various  measures  of  effectiveness 
vary  with  changes  in  those  parameters.  This  will  enable  us  to  determine  the 
parameters  that  we  need  to  estixnate  very  accurately. 

We  have  assumed  that  the  failure  distribution  is  Weibull.  Therefore,  the 
parameters  of  interest  are  the  shape  parameter  a  and  the  scale  parameter  X, 
Besides  these,  other  parameters  are  the  expected  costs  and  the  downtime  of 
preventive  and  failure  replacement.  In  order  to  see  the  changes  in  the 
measures  of  effectiveness  we  hold  all  the  other  variables  constant  and  only 
vary  the  parameter  of  interest.  It  would  be  possible  to  use  Experimental 
Design  techniques  such  as  the  Factorial  Designs  to  study  the  main  effects  on 
the  respective  measures  of  effectiveness.  This  may  help  us  to  identify  the 
parameters  that  are  more  sensitive  and  shall  be  estimated  accurately. 

1.  Alpha 

Table  3.2  shows  the  optimal  value  of  tp*  for  the  minimum  cost  rate, 
the  maximum  availability  and  the  reliability  of  0.95  for  sustaining  a  mission 
of  duration  24  hours  at  age  t;  and  the  respective  optimal  measure  of 
effectiveness  for  various  values  of  the  shape  parameter  alpha.  The  graphical 
plots  are  shown  in  Figures  3.6,  3.7  and  3.8  for  the  cost  rate,  availability  and 
the  reliability  functions  respectively.  It  is  observed  that  as  the  value  of  o  is 
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TABLE  3.2  OPTIMAL  REPLACEMENT  AGE  FOR  VARIOUS 

VALUES  OF  ALPHA 


COST 

AVAiLABniiry 

RELIABILITY 

ALPHA 

V 

$/hr 

t; 

A 

V 

R 

2.6 

1691.8 

29.62 

1228.2 

0.9882 

1647.4 

0.96 

2.8 

1626.0 

29.23 

1166.6 

0.9886 

1424.6 

0.96 

3.0 

1463.6 

28.96 

1126.4 

0.9888 

1371.1 

0.96 

3.2 

1396.3 

28.67 

1104.9 

0.9890 

1331.8 

0.96 

3.6 

1340.7 

28.24 

1083.4 

0.9893 

1289.7 

0.96 

incr«ased  the  replacement  age  geti  chorter  in  order  to  minimize  the  coeti  and 
maximize  the  availability.  This  ii  alao  obvious  from  Figure  3.2  where  the 
hazard  increases  as  alpha  increases  and  as  such  the  failure  rate  also  increases, 
which  in  turns  requires  replacement  action  early  so  as  to  optimize  the 
respective  measures  of  effectiveness.  From  Figure  3.7  it  is  observed  that  the 
curves  are  much  flatter  for  low  values  of  a  and  for  a  ■  3.5  the  curve  falls  quite 
rapidly  on  both  sides  of  the  optimal  value.  Figure  3.8  shows  that  if  we  want  a 
minimum  reliability  of  0.974  then  the  optimal  replacement  age  is 
approximately  980  hours  and  at  this  value  of  reliability  it  is  insensitive  to  the 
shape  parameter  alpha. 
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Tin*  Calf's 


CM  10003 


Figure  8.6  Plot  of  the  Coit  Rate  FunotioB  for  Valuee  of  a. 


CM  0.01) 


TIb»  Chr3 _ 

Figure  8.7  Plot  of  the  Availability  Funotion  for  Valuee  of  a. 


Figure  8.8  Plot  of  the  Reliability  Function  for  Values  of  Alpha. 


From  Table  3.2,  it  is  observed  that  if  the  true  value  of  the  shape 
parameter  a  was  3.2  and  in  the  analysis  we  estimated  it  to  be  3.0,  then  we 
would  have  lost  an  availability  of  0.0002  and  would  have  incurred  additional 
cost  of  $0.28  per  hour.  On  the  other,  if  we  have  estimated  it  to  be  3.6  then  we 
would  have  gained  an  availability  of  0.0003  at  a  lower  cost  by  $0.43  per  hour. 
Therefore  it  is  much  better  to  estimate  the  shape  parameter  higher  than  lower 
and  hence  carry  out  the  replacement  or  overhaul  action  earlier  than  later. 

2.  Lambda 

Table  3.3  shows  the  optimal  value  of  tp*  and  the  respective  optimal 
values  for  the  cost,  availability  and  0.96  reliability  to  sustain  a  mission  of  24 
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hours  duration.  The  graphical  plots  are  shown  in  Figures  3.9, 3.10  and  3.11  for 

the  cost  rate,  availability  and  the  reliability  functions  respectively. 

TABLE  8.8  OPTDf  AL  REPLACEMENT  AGE  FOR  VARIOUS  VALUES 
OF  LAMBDA  (NOTE:  MU  -  1/LAMBDA) 

COST  AVAILABILiry  RELIABILITY 


1360 

1422.2 

29.68 

1102.1 

0.9886 

1326.7 

0.96 

1380 

1443.0 

29.16 

1118.3 

B 

1366.3 

0.96 

1390 

1463.6 

28.95 

1126.4 

0.9888 

1371.1 

0.95 

1400 

1463.9 

28.76 

1134.5 

0.9889 

1386.1 

0.95 

1420 

1484.7 

28.35 

1160.7 

0.9890 

1416.2 

0.96 

From  the  above  table  it  is  observed  that  the  optimal  replacement  age 
for  achieving  a  reliability  of  0.96  when  subjected  to  a  mission  of  duration  of  24 
hours  increases  as  p  increases.  This  is  because  the  higher  the  value  of  p,  the 
reliability  becomes  better.  Also  when  the  value  of  p  increases  the  optimal 
replacement  age  becomes  longer  for  both  the  cost  rate  and  the  availability.  The 
higher  the  value  of  p  the  failure  rate  for  the  Weibull  distribution  decreases  and 
therfore  the  optimal  replacement  age  becomes  longer.  A  small  difference  in  the 
value  of  p  do  not  affect  the  replacement  age  drastically.  However  this  is  not 


Figure  8.9  Plot  of  the  Coet  Kate  Function  for  Values  of  ji. 


Figure  8.10  Plot  of  the  Availability  Function  for  Values  of  |i. 
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Figure  8.11  Plot  of  the  Reliability  Funotloii  for  Valuei  of  )i. 


true  for  the  case  of  alpha.  Therefore,  the  shape  parameter  a  need  to  be 
estimated  more  accurately  than  the  scale  parameter  “K. 

8.  Costs 

The  costs  for  the  failure  replacement  and  the  preventive  replacement 
only  affects  the  cost  rate  function.  These  costs  estimates  are  quite  easily 
available  and  they  are  usually  fairly  accurate  based  on  previous  data.  Since  we 
have  not  discounted  costs  and  cash  flow  problems  in  the  future  it  is  necessary 
to  observe  the  effects  of  the  costs  on  the  cost  rate  function.  It  is  assumed  that 
the  increase  in  costs  are  proportional  such  that  the  ratio  of  replacement  costs, 
C/Cp,  is  the  same.  It  is  observed  that  the  replacement  a^^e  does  not  vary  as 
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long  as  the  ratio  of  the  replacement  costs,  C/Cp,  remains  the  same.  Equation 
(3.7)  dearly  shows  that  if  the  ratio  of  C/Cp,  is  the  same,  there  is  no  effect  on 
the  optimal  replacement  age  tp*.  However  the  cost  per  unit  time  increases 
proportionally  as  the  replacement  costs  increases. 

4.  Downtime 

The  availability  function  greatly  depends  on  the  expected  downtime 
of  the  failure  replacement  (Rf)  and  preventive  replacement  (Rp)  assuming  the 
parameters  for  the  Weibull  distribution  are  estimated  accurately.  In  our 

example  we  have  taken  Rf  a  16  hours  and  the  Rp  a  8  hours.  Because  these 

S’ 

figures  are  small  compared  to  the  optimal  replacement  age  the  availability  is 
very  high.  In  reality  a  failure  at  sea  may  take  days  to  repair  taking  all  the 
logistics  and  administrative  requirements  into  considerations,  espedally  when 
the  ship  is  not  accompanied  by  any  auxilliary  vessel  or  support  ship.  That  is 
the  reason  all  vital  equipments  onboard  have  redimdancies  incorporated.  We 
can  carry  spares  onboard  but  may  not  have  other  support  elements  to  rectify 
the  defect.  Usually  the  downtime  for  a  preventive  replacement  is  constant 
because  it  is  done  at  the  base  with  all  the  support  elements  and  the  work  is 
repetitive  following  standard  procedures.  Therefore,  in  our  analysis  we  shall 
only  vary  R^  to  see  its  effect  on  the  availability.  Table  3.4  shows  the  values  of 
the  availability  if  we  carry  out  the  preventive  replacement  at  1126  hours.  It 
is  observed  that  the  availability  of  the  system  drops  as  R^  increases,  but  is  still 
quite  high. 
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TABLE  8.4  AVAlLABILnY  FOR  VARIOUS  VALUES  OF 

R|.ATTp-1126. 


Expected  Downtime 

Rf  (hours) 

Availability 

A 

24 

0.985 

48 

0.975 

72 

0.965  1 

96 

0.956  1 

C.  EFFECTS  OF  SIMPLE  PREVENTIVE  MAINTENANCE 

In  the  above  policy  it  was  assumed  that  simple  preventive  maintenance 
actions  did  not  improve  the  operational  reliability  of  the  system.  In  real  life  it 
is  evident  that  these  preventive  maintenance  actions  may  definitely  enhance 
the  system  condition  and,  hence  the  reliability,  but  it  will  not  restore  the 
system  to  the  original  state  i.e.  to  be  as  good  as  new.  This  was  also  emphasised 
in  [Ref.  14].  However  in  reality  there  are  situations  for  which  preventive 
maintenance  actions  can  degrade  the  system,  that  is  by  imperfect  repair.  Here 
we  will  concentrate  on  those  actions  that  will  improve  the  condition  by  a 
certain  factor.  So  now  we  shall  incorporate  these  preventive  maintenance 
actions  into  our  model.  The  questions  to  be  asked  is  that  when  should  we  do 
these  type  of  maintenance?  How  do  we  determine  the  improvement  factor? 
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We  can  schedule  the  simple  preventive  maintenance  actions  whenever  the 
operational  reliability  hits  a  minimum  acceptable  limit  say  0.76  and  then 
as  mentioned  above  our  reliability  will  improve  by  a  certain  factor,  thus 
making  the  system  ’younger*  but  not  necessarily  new.  Or  we  could  also  carry 
out  a  planned  replacement  when  the  operational  reliability  hits  the  value 
and  bring  the  condition  of  the  system  back  to  new.  The  appropriate 
maintenance  action  to  be  taken  at  these  critical  times  depends  on  the  measure 
of  effectiveness  required. 

[Ref.  16]  solved  this  problem  by  minimizing  the  expected  cost  rate  as  the 
measure  of  effectiveness.  The  ntunber  of  simple  preventive  maintenance 
actions  before  a  planned  replacement  action  which  gives  a  minimum  cost  rate 
is  then  obtained  by  evaluating  the  cost  rate  whenever  the  system  reaches  R^ua. 
Here  we  shall  carry  out  the  appropriate  maintenance  actions  so  as  to  maximize 
the  availability  of  the  system. 

Notation  is  as  follows: 

IP  simple  preventive  maintenance 

2P  planned  preventive  repla(»ment 

2C  failure  replacement 

T  time  to  failure 

Ri(t))  probability  of  no  failure  during  (t|.i,t|) 

Fi(t|)  probability  of  a  failure  during  (t|.i,t,) 

F(t)  probability  distribution  hmction  of  the  time  to  failure 
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f(t)  probability  density  function  of  the  time  to  failure 
h(t)  failure  rate  function  of  the  time  to  failure 
Ry  time  required  for  each  maintenance  action  iMl,2  u«P,C 
Cy  cost  of  each  maintenance  action  ial,2  ‘J«P,C 
Rnia  minimum  acceptable  operational  reliability  limit 
t|  ith  time  the  system  reaches 

Let  I  be  the  improvement  factor,  then  if  maintenance  type  IP  is  carried 
out  at  time  t^  then  the  system  age  is  reduced  firom  tj,  to  ( t^I )  and  therefore 

t,  -  ti  +  ti  (1  -  1) 

and  t,  -  t,  +  <t,  -  ti)  (I  -  1) 

in  general  t.  -  t,.i  +  (1  -  i)“-‘  tj  (8.18) 

As  I  -4  1,  the  effect  on  the  system  approaches  "bad  as  old"  and 
As  I  0,  the  effect  approaches  "good  as  new"  [Bef.  16]. 

The  value  of  t(  can  be  computed  as  follows: 


r,  -  (3.19) 

and  firom  equation  (3.19)  we  can  compute  t^  for  i  ■>  2,3,...,n. 

The  value  of  I  can  be  estimated  from  past  records  from  the  data  collected 
on  similar  equipment  where  performance  measurements  are  taken  before  and 
after  the  simple  preventive  maintenance  actions.  In  some  cases  performance 
measurement  techniques  such  as  condition  monitoring  by  vibration  analysis 
are  employed.  This  is  largely  used  in  rotating  machinery.  As  mentioned  earlier, 


40 


simple  preventive  maintenance  actions  constitutes  cleaning,  adjustments, 
replacing  worn  out  parts,  tuning  and  minor  calibration.  All  these  actions 
normally  improve  the  reliability  which  is  a  fimction  of  the  failure  rate.  So 
these  actions  xoay  prolong  the  life  of  the  system,  and  enable  the  system  to  be 
'younger*.  From  experience,  it  is  found  that  as  the  system  becomes  older  these 
simple  preventive  maintenance  actions  becomes  ineffective  and  as  such  the 
system  may  require  a  complete  overhaul  which  on  completion  is  as  good  as 
new.  For  our  case  here  we  shall  consider  that  on  completion  of  each  simple 
preventive  maintenance  action,  the  system  age  is  reduced  from  t^  to  t^  •  (t|  •  t|.i) 
(1  - 1)  where  ti  is  the  time  of  the  simple  preventive  maintenance  action. 

A  cycle  is  completed  each  time  a  replacement  takes  place  which  could  be 
either  due  to  a  failure  replacement  or  a  planned  replacement.  On  completion 
of  the  replacement  action  the  cycle  probabilistically  starts  all  over. 

1.  Case  1 

T  >  ta  (Failure  occurs  after  the  type  2P  maintenance) 

IP  IP  2P 

I - ^ ^ - 1 

0  ti  t,  . .  .  t. 

The  costs  associated  are: 

a)  Expected  costs  of  type  IP  maintenance. 

If  planned  replacement  takM  place  at  to,  then  there  will  be  (n  •  1) 
maint  nanco  actions  of  type  IP  and  each  costs  C^p. 
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Then  total  cost  s  (n  -  1)  Cxp 

But  this  event  will  occur  if  there  are  no  failures  in  the  first  (n-1) 
intervals.  The  Probability  of  this  event  is  ■  Ri(ti)Ra(ta)..,Ro(tB).  It  is  assumed 
that  the  failures  are  statistically  independent.  The  failure  in  each  interval  is 
independent  from  the  failure  in  next  interval. 

Therefore  expected  costs  of  type  IP  maintenance  ■ 

in  -  !)•  n  m  (3*20) 

where  R|(t|) »  1  •  Probability  of  a  failure  in  the  (n*l)  intervals  and  is  given  by 

m  ■  1  “  TO  -  “  Vi))]  (3-21) 

b)  Expected  costs  of  type  2P  maintenance. 

Since  there  will  he  only  one  Type  2P  maintenance  in  an  interval,  the 
expected  costs  of  type  2P  maintenance  ■* 

Car'flW  (3-22) 

1-1 

The  expected  downtime  can  be  obtained  by  substituting  the  cost  by 
downtime  in  equations  (3.19)  and  (3.20). 

The  downtime  associated  are: 

a)  Expected  downtime  of  type  IP  maintenance  ■ 

in  -  (3.23) 

<-i 

b)  Expected  downtime  of  type  2P  maintenance  ■ 
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(3.24) 


^•n  w 

2.  CsM  2 

t.-i  <  T  <  t. 


» 


0 


IP 

■f- 

tx 


IP 


^#“1 


4^ 


The  costs  associated  are: 

a.  Expected  cost  of  type  IP  maintenance 

If  the  flsilure  occurs  between  t^,,  and  then  there  will  be  (n-1) 
preventive  maintenance  actions  and  the  Probability  of  this  event  is 
■  Probability  of  no  failures  in  the  first  (n>l)  intervals  multiplied  by  the 
Probability  of  a  failure  between  t^fi  and  t^.  Again  we  assume  that  failures  are 
independent.  Then  tlie  expected  cost  of  type  IP  maintenance  ■ 


l-l  >-l 


(3.25) 


b.  Expected  costs  of  type  2C  maintenance  ■ 

El^ac-iiwWl  (3'“) 

M  /-I 

The  downtime  associated  are: 

a.  Expected  downtime  of  type  IP  maintenance  ■ 
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(347) 


EW“ 

1-1  ji-i 

b.  Expected  downtime  of  type  2C  maintenance  ■ 

g[^*n  W  ‘W  C3.28) 

The  expected  operation  timo  until  t^  is  the  Bum  of  the  average 
operation  time  until  t^.,^  and  the  average  operation  time  during  (to.i»to). 

The  expected  downtime  is  the  sum  of  all  the  downtimes  in  equations 
(3.23),  (3.24),  (3.27)  and  (3.28). 

From  our  earlier  discussion  we  know  that  the  mean  life  of  the 
system,  assuming  only  1  interval,  is  given  by 

«(»,)  ■  Q  H  -  m  *  osm 

and  the  above  equation  can  also  be  written  as 

M(t.)  -  dt  *  (3.30) 

and  for  n  intervals  the  equation  can  be  written  as 

slnw - fw)i| 

1-1  I  ,1-1  *  Vi'S!*! 
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Total  Oparatlon  Tima 

Mow  Availability  - - 

Total  Qparation  Tima  -i-  Total  Downtima 

It  vrill  be  shown  later  that  if  the  expected  downtime  for  a  type  IP 
maintenance  (simple  preventive  maintenance)  is  small  then  the  availability 
function  given  by  the  above  equation  is  a  good  approximation. 

The  expected  costs  is  the  sum  of  all  the  costs  in  equations  (3.20), 
(3.22),  (3.26)  and  (3.26).  Then  the  cost  rate  function,  C(tp),  is  given  by 

Total  Ixpaotad  Costa 

C(tp)  - - 

Total  Oparation  Tima  +  Total  Downtiaw 
8t  Numaiioml  Illustration 

Now  let  us  take  a  case  example  with  the  following  data: 

Let  the  time  to  failure,  T,  follow  a  Weibull  distribution  with  shape 
parameter  a  and  scale  parameter  X. 

Lot 

a  ■  3,0 
X  m  1/1390 
Rip  ■  1  hour 
Rap  ■  8  hours 
Rao  48  hours 

Using  the  above  data  and  equations,  a  Fortran  program  was  written 
[Annex  A]  to  determine  the  optimum  times  to  carry  out  simple  preventive 
maintenance  actions  cmd  complete  overhauls  or  preventive  replacement  so  as 
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to  maximize  the  availability.  The  results  are  tabulated  in  Tables  3.6,  3.6,  and 
3.7  for  values  of  I  a  0.1,  0.6  and  0.9  respectively. 

When  I  a  0.1,  carry  out  simple  preventive  maintenance  actions  at 

times: 

•  917.6  hours 

•  1743.4  hours 

•  2486.7  hours 

and  preventive  replacement  or  complete  overhaul  at  time  3166.6  hours  to 
achieve  a  maximum  availability  of  0.9836. 

When  I  a  0.6,  carry  out  simple  preventive  maintenance  action  at  time 
917.6  hours  and  a  preventive  replacement  at  time  1376.4  hours  to  achieve  a 
maximum  availability  of  0.9802. 

When  I  a  0.9.  do  not  carry  out  any  simple  preventive  maintenance 
action  and  at  time  917.6  hours  replace  or  carry  a  complete  overhaul  of  the 
system  so  as  to  achieve  a  maximum  availability  of  0.9794. 

It  is  observed  that  as  the  improvement  factor,  I  -4  0,  the  system 
becomes  "good  as  new"  after  each  simple  preventive  maintenance  action  and 
as  the  improvement  factor,!  1,  the  system  becomes  "bad  as  old"  after  each 
preventive  maintenance  action. 
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TABLE  8.5  REPLACEMENT  TIMES  FOR  IMPROVEMENT  FACTOR, 
I  >  0.1  AND  Rv  -  8  HOURS,  Rm  -  48  HOURS,  Rip  -  1  HOUR,  o  -  3.0, 

|i  -  1800  HOURS. 


N 

TIME  (HOURS) 

AVAILABILITY 

1 

917.6 

0.979421 

2 

1743.4 

0.982197 

3 

2486.7 

0.983306 

4 

3166.6 

0.983691 

6 

3767.7 

0.983686 

6 

4299.6 

0.983480 

7 

4787.1 

0.983344 

8 

6226.0 

0.983206 

9 

6621.0 

0.983076 

10 

6976.6 

0.982966 

11 

6296.6 

0.982848 

12 

6684.4 

0.982761 

13 

6843.6 

0.982664 

14 

7076.8 

0.982687 

16 

7286.7 

0.982618 
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TABLE  3.6  REPLACEMENT  TIMES  FOR  IMPROVEMENT  FACTOR, 
I  -  0.5  AND  R«.  -  8  HOURS,  R^  -  48  HOURS,  R^p  -  1  HOUR,  a  -  8.0, 

H  m  1390  HOURS. 


TIME  (HOURS) 

AVAILABILITY 

917.6 

0.979421 

1376.4 

0.980166 

1606.8 

0.980074 

1720.6 

0.979848 

1777.8 

0.979644 

1806.6 

0.979483 

1820.9 

0.979368 

1828.0 

0.979261 

1831.6 

0.979184 

1833.4 

0.979122 

1834.3 

0.979071 

1834.7 

0.979029 

1836.0 

0.978994 

1836.1 

0.978964 

16 


1836.1 


0.978938 


TABLE  3.7  REPLACEMENT  TIMES  FOB  IMPROVEMENT  FACTOR, 
1  -  0.&'  AND  R«.  n  6  HOURS,  Hk  ■  48  HOURS,  Rip  ■  1  HOUR,  a  ■  3.0, 

11-1890  HOURS. 


N 

TIME  (HOURS) 

AVAILABH^IIY 

1 

917.6 

0.979421 

2 

1009.4 

0.979042 

3 

1018.6 

0.978466 

4 

1019.6 

0.977966 

6 

1019.6 

0.977640 

6 

1019.6 

0.977198 

7 

1019.6 

0.976899 

8 

1019.6 

0.976646 

9 

1019.6 

0.976426 

10 

1019.6 

0.976123 
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IV.  POUCY II  (MINIMAL  REPAIR  WITH  AGE  REPLACEMENT) 

As  an  alternative  assume  that  the  system  consists  of  several  components 
and  the  system  fails  when  one  of  the  components  fail,  that  is  components  are 
connected  in  series.  We  shall  assume  tide  components  have  independent  IFR 
distribution.  Then  the  system  lifetime  is  also  IFR.  (The  IFR  property  is 
maintained  under  the  formation  of  the  series  system).  Preventive  replacement 
is  carried  out  when  the  system  reaches  the  age  t^.  Between  the  preventive 
maintenance,  failures  are  repaired  as  quickly  as  possible  (minimal  repairs) 
either  by  replacing  the  failed  component  with  a  new  one  or  repairing  the  failed 
part.  We  also  assume  that  the  system  failure  rate  is  not  disturbed  on 
completion  of  the  minimal  repair  i.e.  the  failure  rate  is  the  same  as  before 
failure.  Tf  the  failure  occurs  at  time  t  <  t^  the  failure  rate  of  the  system  just 
after  the  minimal  repair  is  h(t).  Failures  are  detected  immediately. 

This  idea  of  minimal  repair  was  first  introduced  by  Barlow  and  it  is 
described  in  [Ref.  7].  After  that,  many  authors  expanded  on  these  ideas  and 
formulated  various  models  associated  with  minimal  repairs;  they  used  the 
expected  cost  rate  as  the  measiure  of  effectiveness.  In  this  policy  we  will 
emphasize  the  availability  of  the  system  as  our  measure  of  effectiveness. 
However  we  will  also  formulate  the  expected  cost  rate  models. 
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Let  0  <:  t  be  the  system  operating  time  since  last  replacement.  Then  the 
probability  of  a  failure  occurring  in  [t,t4>dt]  is  hCtIdt  where  h(t)  is  the  failure 
rate  of  the  system. 

Let  N(tp)  ■  number  of  failures  occuring  in  time  (0,tp) 

When  repairs  are  minimal,  { N(tp) ,  tp  ^  0  }  is  a  Nonhomogeneous  Poisson 
Process  with  intensity  iUnction  h(tp). 

The  probability  that  the  system  will  experience  n  failures  in  the  interval 
(0,tp)  is  given  by 

(4..) 

where 


m 


(V  -  «i  MV  1  ■ 


m(tp)  is  called  the  mean  value  function. 


(42) 


A.  COST  MODEL 

The  age  t^*  at  wliich  the  operating  system  is  replaced  depends  on  the 
following  factors: 

•  Failure  distribution  (IFR). 

•  Costs  of  minimal  repair  and  preventive  replacement. 

•  Downtime  of  minimal  repair  and  preventive  replacement. 
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*  Measure  of  effectiveness; 

-  Minimize  costs. 

-  Maximize  availability. 


We  have  assumed  that  there  are  no  simple  preventive  maintenance 
actions. 


-Re- 


* 


failure 


Let  Cf  be  the  mean  cost  of  each  minimal  repair. 

Cp  be  the  mean  cost  of  preventive  replacement. 

Rp  be  the  mean  time  of  preventive  replacement. 

Rr  be  the  mean  time  of  a  minimal  repair, 
and  tp*  be  the  planned  replacement  age. 

Cf  >  Cp  and  Rf  >  R^ 

Then,  assuming  that  the  cost  of  downtime  is  negligible,  the  total  cost  is 
given  by 


Mg 

C  -  E  C;.  *  (4J) 

where  Cf,  is  the  cost  of  the  i*th  minimal  repair  and  the  length  of  the  cycle  is 
given  by 

L  «  tp  +  K,  (4.4) 
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By  using  equation  (3.1)  and  assuming  that  the  time  for  a  minimal  repair 
is  very  small  compared  to  that  of  the  length  of  a  cycle,  the  expecte'^.  cost  rate 
function  is  given  by 


C(rp 


c,  *»  *  ♦  c. 


(4.5) 


and  we  wish  to  minimize  the  cost  per  unit  time,  so  we  set  the  derivative  of 
equation  (4.6)  to  0  and  we  obtain 

C. 


wpii,  *  jy  -  /^'  *(#  d(  - 


(4« 


Then  the  value  of  tp*  that  satisfies  the  above  equation  is  the  optimal 
replacement  age. 

For  a  time  to  failure,  T,  following  a  Weibull  distribution  with  shape 
parameter  a  and  scale  parameter  X  the  expected  cost  rate  function  is  given  as 
below 


C(y  .  (4.7) 

For  given  values  of  Cf ,  Cp ,  R, ,  Rf  and  the  parameters  a  and  X  we  can 
find  the  optimal  values  of  tp*  that  minimize  the  expected  cost  rate. 

Now  let  us  take  a  case  example  with  the  data  as  follows,  but  first  without 


preventive  maintenance  action: 

Cp  -  $  26000 
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Cf-$  1000 
Rp  ■  8  hours 
Rf  ■  1  hour 
a  ■  3.0 

X  ■  1/1390  hours 


Figure  4.1  is  a  plot  of  the  cost  rate  function  as  in  equation  (4.7)  and  the 
optimal  value  of  tp*  that  minimizes  the  cost  rate  fiuiction  is  3222  hours  at  a 
costof  $11.60/hour. 


Table  4.1  shows  the  values  of  the  cost  rate  for  the  replacement  age  t^ 
close  to  the  optimal  value.  From  Figure  4.1  it  is  observed  that  the  curve  is 
fairly  flat  near  the  optimal  value  thus  giving  some  flexibility  to  the  decision 
maker. 
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TABLE  4.1  VALUES  OF  COST/HOUR  FOR  VARIOUS 
REPLACEMENT  AGE. 


Replacement  Age 

Coat 

(houra) 

(4/hour) 

11.66 

3060.0 

11.63 

3100.0 

11.61 

3160.0 

11.60 

3200.C 

11.60 

3260.0 

11.60 

3300.0 

11.60 

3360.0 

11.61 

3400.0 

. 

11.63 

3460.0 

11.66 

3600.0 

11.68 

Now  suppose  there  are  s  identical  components  which  operate  and  fail 
independently  and  have  the  same  failure  distributions  with  the  same 
parameters.  All  of  these  components  are  replaced  at  time  tp  at  a  cost  of  Cp . 
When  each  of  the  component  fails  it  ie  replaced  individually  or  undergoes 
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minimal  repair  without  affecting  the  system  failure  rate  as  a  whole.  Then  the 
expected  cost  rate  ffinction  becomes 


cap 


*c. 


(4.8) 


Now  let  us  imagine  that  we  cany  out  a  simple  preventive  maintenance 
action  at  some  time  tp  and  after  N  of  these  preventive  maintenance  actions  are 
carried  out  we  either  replace  the  system  or  carry  out  a  complete  overhaul 
which  on  completion  is  like  new.  Any  failures  between  the  preventive 
replacement  are  treated  as  minimal  failures  and  are  repaired  quickly  in 
negligible  downtime  at  a  cost  of  Cf.  The  cost  of  each  simple  preventive 
maintenance  action  is  0,p.  It  is  assumed  that  alter  each  simple  preventive 
maintenance  action  the  system  improves  by  a  certain  factor,  so  for  simplicity 
we  say  that  the  system  becomes  ’younger*.  Nakagawa  [Ref.  16].  This  means 
that  if  t  is  the  time  of  a  preventive  maintenance  action  the  failure  rate  on 
completion  is  h(t-x)  where  x  is  the  amount  of  time  by  which  the  system  has 
become  ’younger*.  However  the  failure  rate  after  a  minimal  repair  stays  the 
same.  The  value  of  x  can  perhaps  be  determined  from  past  historical  records 
usixig  performance  measurements  techniques  such  as  condition  monitoring  by 
vibration  analysis  or  some  output  parameter.  However  a  methodology  for 
characterizing  the  effective  age  reduction,  x,  remains  to  be  developed. 

According  to  the  above  assumption  the  replacement  age  is  Ntp,  and 
there  are  N*1  simple  preventive  maintenance  actions. 


S6 


A  cycle  is  completed  at  (Nt^  4*  11^)  aiid  the  process  probabilistically 
starts  all  over. 


0 


I 


f  — -  minimal  failure 

Rp  --  expected  time  of  preventive  replacement 

t),  i  m  simple  preventive  maintenance 

Ntp  preventive  replacement 

(  Note:  At  ti  the  age  of  the  system  is  (ti  -  x) ) 

The  total  expected  costs,  E[C],  incurred  in  a  cycle  is  given  by 

-  C,  +  C^+  C,  4  +  ...  4  C,  (4.9) 

for  0  ^  X  :S  tp 

and  this  equation  can  be  simplified  to 

*IC!  -  C,  £  A<0  *  *  (AT-DC,  *  C,  (4.1(0 

for  0  ^  X  $  tp 

and  then  the  expected  cost  rate  fUnctionc  C(tp),  is  given  by 

•  CCl 

C(tp)  - - 

l[L] 


where  E[L]  ■  Ntp  +  Rp. 
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For  a  time  to  failure,  T,  following  a  Weibull  distribution  with  shape 
parameter  a  and  scale  parameter  X  he  expected  cost  rate  hmcUon  is  given  as 
follows 

9  E  <»  ■*  *  c, 

- - -  (4.11, 

When  C,p  ■  0  and  x  ■  0,  that  is  there  is  no  preventive  maintenance 
action,  then  the  above  equation  is  the  same  as  equation  (4.7). 

Now  we  shall  plan  to  carry  out  a  preventive  maintenance  at  some 
time  tp  at  a  cost  of  C,p  ■  $  600.00.  On  completion  of  this  simple  preventive 
maintenance  action,  the  age  of  the  system  becomes  younger  by  some  value  x. 
When  X  ■  0,  it  is  observed  that  there  is  no  improvement  to  the  optimal  planned 
replacement  age  but  we  have  incurred  an  additional  cost  at  $11.76/hour 
instead  of  the  original  cost  at  $1 1.60/hour.  As  x  increases  firom  0,  the  optimal 
planned  replacement  age  increases  and  also  the  cost  rate  reduces.  This  is  for 
the  case  of  carrying  out  only  one  preventive  maintenance  in  between  the 
planned  replacement  or  overhaul  age.  We  could  carry  on  the  same  analysis 
incorporating  more  preventive  maintenance  actions  and  the  results  will  give 
more  improvement.  However  this  greatly  depends  on  the  parameters  such  as 
the  value  of  x  and  the  cost  of  preventive  maintenance  assuming  the  other 
parameters  remain  fixed. 
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In  the  real  scenario  as  the  system  gets  older  there  is  more 
deterioration  and  on  completion  of  each  preventive  maintenance  action  the 
value  of  X  probably  reduces  over  time;  that  is;  it  is  some  function  of  t,  and 
there  comes  a  stage  at  which  preventive  maintenance  action  will  not  improve 
the  system  any  longer.  For  simplicity,  however  we  assume  that  x  is  constant. 
It  should  be  noted  that  there  will  be  values  of  age  for  which  (age  •  x)  <  0  ,  we 
therefore  take  the  maximum  (0.  (age  -  x))  as  the  age  uf  the  system  on 
completion  of  a  simple  preventive  maintenanoe  action. 

The  Figure  4.2  and  Table  4.2  illustrate  that  simple  preventive 
maintenance  actions  increases  the  optimal  replacement  age. 

Table  4.2  shows  the  time  for  a  preventive  replacement  (TPR)  and  the 
time  for  a  simple  preventive  maintenanoe  (TPM)  for  various  values  of  x. 

B.  AVAILABILITY  MODEL 

Now  let  us  take  the  availability  as  our  measure  of  effectiveness.  From 
earlier  discussion  Availability  is  given  by 

Mean  Lif<s  During  A  Cycle 

A(t^)  - - 

Vetal  length  of  cycle 


X.)d. 


(4.12) 


and  we  wish  to  maximize  the  availability  function,  so  we  set  the  derivative  of 
equation  (4.12)  to  0  and  we  obtain 
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Figure  4.2  Plot  of  the  Coat  Funotioii  with  Preventive  Replacement 
for  Varioua  Valuea  of  x.  Cgp  ■  $500.00. 


TABLE  4.2  OPTIMAL  PREVENTIVE  REPLACEMENT  AGE  WITH 
PREVENTIVE  MAINTENANCE  FOR  VARIOUS  VALUES  OF  X. 

Cgp  a  $600.00. 


z 

(hours) 

TPR 

(hours  a  1) 

TPM 

(hours  ±  0.5) 

Cost 

($/hour) 

0 

3244 

1622 

11.76 

100 

3281 

1641 

11.48 

200 

3320 

1660 

11.22 

600 

3443 

1722 

10.49 

800 

3673 

1787 

9.83 

1000 

3667 

1834 

9.43 

(4.13) 


Then  the  value  of  tp*  that  Badsfiee  the  above  equation  is  the  optimal 
replacement  age. 

For  a  time  to  failure,  T,  following  a  Weibull  distribution  with  shape 
parameter  a  and  scale  parameter  X  the  availability  function  is  given  as  follows 

The  availability  ftuicdon  in  equation  (4.12)  is  an  approximation  because 
we  do  not  take  into  account  of  the  downtime  during  a  minimal  repair,  which 
in  reality  is  not  entirely  valid.  In  the  equation  we  assume  that  the  system  can 
still  fail  when  the  system  is  down  (we  are  integrating  over  (0,tp)  and  some  of 
that  is  downtime).  However,  this  approximation  is  quite  accurate  if  the 
downtimes  are  small,  as  is  likely  to  be  true  in  practice. 

We  can  find  an  exact  solution  to  equation  (4.12)  as  follows: 

Let 

A(t)  m  Probability  the  system  is  available  at  time  t  (age  t) 
following  the  last  replacement. 
h(t)  ■  failure  rate  at  time  t. 
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p(t)  =  ropair  rate  at  time  t.  For  our  case  this  is  assiimed  to  be 
constant  and  denoted  by  p.  Then  p  s  It  ^ 

(It  should  be  noted  that  it  may  take  longer  to  repair  an  older  system) 
Then, 

A<t+4t)  -  A(b)(l  -  2i(t)dt)  +  {1  -  A{t))p(t)dt  +  o(dt)  (4.15) 
where  o(dt)  represents  higher  order  terms  which  are  s  0 
From  equation  (4.16),  arranging  the  terms  we  get 


A(t+dt)  -  A(t) 


-  -h(t)  A(t)  +  p{t)  (1  -  A(t)) 
dt 

(4.16) 

dA(t) 

-  -h(t)  A(t)  -  p(t)  A(t)  +  p(t) 
dt 

(4.17) 

dA(t) 

+  (h(t)  +  p(t))  A(t)  •  p(t) 
dt 

(4,18) 

^  M(f)  ^  -  m  ^ 

dt 

(4.19) 

,  m .  m ^ 

(4.20) 

(4.21) 

For  a  time  to  failure,  T,  following  a  WeibuU  Distribution  with  shape  parameter 
a  and  scale  parameter  X,  the  failure  rate  h(t)  is  given  by 

h(t)  •  «>l,l>0/i0  (4.22) 

then  subtituting  A(0}  s  1  and  ,  p(t) «  (Rf)'^ ;  A(t)  can  be  simplified  as  follows: 
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A(r)  -  € 


1  -K4«)“  -  ai,)*  - -^1 

’^io*  '*■ 


t  .  h. 


But  the  integral 


I  •  ^ 


has  no  closed  form 


solution.  Therefore  we  have  to  evaluate  it  nmnerically. 

Using  the  same  principle  as  that  in  [Ref.  17]  the  average  availability  over 
a  (^cle  of  length  tp  Rp  and  hence  in  the  long  run  is  given  by 

■  rfr  /o’ 

By  taking  the  derivative  of  equation  (4.19)»  we  obtain 


d4(fp)  '■  ^  A(t)  dt 

A,  ‘ 


(4.25) 


Setting 

dX(tp) 

-  ■  0  ,  for  tp  >  0  then  tp  is  a  candidate  for 

d(tp)  a  time  between  the  end  of  one  preventive 

replacement  and  the  beginning  of  the  next 

and  tp*  is  the  value  that  satisfies  the  following  equation 


Ad,-)  -  — i — /y  Ad)  It  (4M) 

To  solve  equation  (4.24)  the  IMSL  subroutine  DQDAG  was  used  to 
compute  the  integral  in  equation  (4.23)  for  various  values  of  t.  So  we  obtained 
A(t)  for  various  values  of  t  and  the  graph  of  this  function  is  shown  in  Figure 
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4.3.  Then  we  used  the  IMSL  subroutine  DCSINT  [Ref.  13]  to  compute  the  cubic 
spline  interpolant  to  the  set  of  data  points  obtained  earlier  (values  of  t  and  the 
corresponding  values  of  A(t) ).  Finally  we  used  the  IMSL  subroutine  DCSITG 
[Ref.  13]  to  evaluate  the  integral  of  the  cubic  spline  for  various  values  of  tp.  A 
brief  description  of  the  IMSL  subroutines  DCSINT  and  DCSITG  are  in 
Appendix  D  and  E  respectively.  The  values  obtained  by  using  this  ’exact’ 
method  wao  compared  with  those  values  obtained  using  the  approximation  as 
in  equation  (4.12).  These  values  are  tabulated  below.  It  is  observed  that  the 


Figure  4.8  Plot  of  the  Availability  Function  for  Valuea  of  aa  8.0, 
|i  ■  1890,  ■  8  hours,  ■  8  hours. 
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both  the  approximate  and  exact  solutions  are  the  same  to  four  decimal  places 
and  the  approximation  gives  an  accurate  solution.  This  is  true  for  small  values 
of  Rf  and  as  Rf  is  increased  from  1.0  hour  to  8  hours,  both  the  approximation 
and  the  exact  solutions  are  the  same  to  two  decimal  places,  as  shown  in  Table 
4.4.  However  the  optimal  replacement  age  are  the  same  for  both  the  cases.  As 
Rf  increases  the  accuracy  of  the  approximation  diminishes. 

Table  4.3  shows  the  values  of  availability  for  a  ■  3.0,  p  ■  1390,  Rf  ■  17 
hours.  Rp  ■  8  hours.  The  optimal  replacement  age  for  the  two  methods  and  the 
maximum  availability  are: 

•  approximate  solution  2203  hours  0.994681 

•  exact  solution  2208  hours  0.994689 

It  is  observed  that  the  approximation  gives  very  accureta  results. 

Table  4.4  shows  the  values  of  availability  for  a  ■  3.0.  p  ■  1390.  Rf  u  8 
hours.  R^  ■  8  hours.  The  optimsil  replacement  age  for  the  two  methods  and  the 
maximum  availability  are: 

•  approximate  solution  1099  hours  0.989201 

•  exact  solution  1108  hours  0.989301 

It  is  again  observed  that  the  approximation  gives  very  accurate  results. 
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Table  4.5  shows  the  values  of  availability  for  a  =  2.6,  p  =  1390,  Rf  s  1 
hour,  Rp «  8  hours.  The  optimal  replacement  age  for  the  two  methods  and  the 
maximum  availability  are: 

*  approximate  solution  2711  hours  0.996104 

•  exact  solution  2717  hours  0.996112 

Once  again  it  is  observed  that  the  approximation  gives  very  accurate 
results.  If  the  actual  shape  parameter  was  2.6  and  if  we  have  estimated  it  to 
be  3.0,  then  we  would  liave  lost  an  availability  of  0.000623  which  can  be 
considered  negligible. 

Table  4.6  shows  the  values  of  availability  for  a  ■  3.6,  p  ■  1390,  Rf  m  1 
hour,  Rp  ■  8  hours.  The  optimal  replacement  age  for  the  two  methods  and  the 
maximum  availability  are: 

•  approximate  solution  1936  hours  0.994243 

*  exact  solution  1940  hours  0.994263 

Once  again  it  is  observed  that  the  approximation  gives  very  accurate 
results.  If  the  actual  shape  parameter  was  2.6  and  if  we  have  estimated  it  to 
be  8.6,  then  we  would  have  lost  an  availability  of  0.000869  which  again  can  be 
considered  negligible.  This  shows  that  the  shape  parameter  does  not  effect  the 
availability  drastically  and  a  close  estimate  is  sufficient. 


66 


Table  4.7  shows  the  values  of  availability  for  a  =  3.0,  =  1360,  Rf  =  1 

hoiu*,  Rp  B  8  hours.  The  optimal  replacement  age  for  the  two  methods  and  the 
maximum  availability  are: 

*  approximate  solution  2140  hours  0.994421 

*  exact  solution  2145  hours  0.994430 

It  is  observed  that  the  approximation  gives  very  accurate  results.  If  the 
actual  scale  parameter  was  1390  and  if  we  have  estimated  it  to  be  1350,  then 
we  would  have  lost  an  availability  of  0.000159  which  again  can  be  considered 
negligible. 

Table  4.8  shows  the  values  of  availability  for  a  ■  3.0,  p  ■  1450,  Rf  ■  1 
hour,  Rp  ■  8  hours.  The  optimal  replacement  age  for  the  two  methods  and  the 
maximiun  availability  are; 

*  approximate  solution  2298  hottrs  0.994805 

*  exact  solution  2302  hours  0.994812 

It  is  observed  that  the  approximation  gives  very  accurate  results.  If  the 
actual  scale  parameter  was  1450  and  if  we  have  estimated  it  to  be  1350,  then 
we  would  have  lost  an  availability  of  0.000382  which  again  can  be  considered 
negligible.  This  shows  that  the  scale  parameter  does  not  effect  the  availability 
drastically  and  a  close  estimate  is  sufficient. 
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TABLE  4.8  COMPARISON  OF  AVAILABILlTy  FOR  VALUES  OF 
a  ■  8.0»  |i  -  1890,  R,  -  1  HOUR,  Rp  -  8  HOURS. 


Replaoement 
Age  (hours) 

Approximate 

Availability 

Exact 

Availability 

1200.0 

0.992846 

0.992847 

1300.0 

0.993268 

0.993261 

1400.0 

0.993693 

0.993696 

1600.0 

0.993862 

0.993866 

1600.0 

0.994076 

0.994080 

1700.0 

0.994245 

0.994249 

1800.0 

0.994374 

0.994379 

1900.0 

0.994469 

0.994474 

0.994532 

0.994639 

2100.0 

0.994669 

0.994676 

2200.0 

0.994681 

0.994689 

2300.0 

0.994671 

0.994680 

2400.0 

0.994540 

0.994661 

2600.0 

0.994490 

0.994603 

2600.0 

0.994423 

0.994437 

TABLE  4.4  COMPARISON  OF  AVAILABILITY  FOR  VALUES  OF 
a  ■  8.0,  |i  ■  1800,  Rr  -  8  HOURS,  Rp  -  8  HOURS. 


Replacement 
Age  (hours) 

Approximate 

Availability 

Exact 

Availability 

800.0 

0.988211 

0.988273 

900.0 

0.988798 

0.988870 

0.989108 

0.989193 

1100.0 

0.989201 

0.989300 

1200.0 

0.989116 

0.989232 

1300.0 

0.988880 

0.989014 

1400.0 

0.988613 

0.988669 

1600.0 

0.988028 

0.988210 

1600.0 

0.987437 

0.987648 

1700.0 

0.986748 

0.986993 

1800.0 

0.986967 

0.986262 

0.984148 

0.984630 

2200.0 

0.982012 

0.982619 

2400.0 

0.979677 

0.980244 

2600.0 

0.976867 

0.977724 

TABLE  4.6  COMPARISON  OF  AVAILABILITY  FOR  VALUES  OF 
a  M  2.6,  |i  -  1890,  Rr  -  1  HOUR,  Rp  -  8  HOURS. 


Raplaoement 
Age  (hours) 

Approximate 

Availability 

Exact 

Availability 

n^QQ^iiiiiiiii 

0.994779 

0.994783 

2100.0 

0.994874 

0.994878 

2200.0 

0.994949 

0.994964 

2300.0 

0.996008 

0.996013 

2400.0 

0.996051 

0.996057 

2600.0 

0.996080 

0.996087 

2600.0 

0.995098 

0.996106 

2700.0 

0.996104 

0,996112 

2800.0 

0.995100 

0.996108 

2900.0 

0.995087 

0.996096 

IIIIIIIQQIjQIIIIII^ 

0.996066 

0.996076 

3100.0 

0.995036 

0.995047 

3200.0 

0.995000 

0.996011 

3300.0 

0.994956 

0.994969 

3400.0 

0.994907 

0.994921 
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TABLE  4.6  COMPARISON  OF  AVAILABILITY  FOR  VALUES  OF 
a  -  8.5,  |i  -  1890,  R^  -  1  HOUR,  Rp  -  8  HOURS. 


Replacement 
Age  (hours) 

Approximate 

AvallabiUty 

Exact 

Availability 

1200.0 

0.992883 

0.992886 

1300.0 

0.993279 

0.993281 

1400.0 

0.993690 

0.993693 

1600.0 

0.993829 

0.993833 

1600.0 

0.994007 

0.994012 

1700.0 

0.994132 

0.994137 

1800.0 

0.994208 

0.994216 

1900.0 

0.994242 

0.994250 

2000.0 

0.994236 

0.994246 

2100.0 

0.994194 

0.994206 

2200.0 

0.994118 

0.994132 

2300.0 

0.994009 

0.994026 

2400.0 

0.993869 

0.993889 

2600.0 

0.993699 

0.993723 

2600.0 

0.993601 

0.993629 

TABLE  4.7  COMPARISON  OF  AVAILABIUTY  FOR  VALUES  OF 
a  -  8.0,  fi  -  1850,  Rr  -  1  HOUR,  Rp  -  8  HOURS. 


Replacement 
Age  (hours) 

Approximate 

Availability 

Exact 

Availability 

1300.0 

0.993201 

0.993204 

1400.0 

0.993626 

0.993629 

1600.0 

0.993786 

0.993789 

1600.0 

0.993990 

0.993993 

1700.0 

0.994147 

0.994162 

1800.0 

0.994264 

0.994269 

1900.0 

0.994346 

0.994362 

0.994397 

0.994404 

2100.0 

0.994419 

0.994428 

2200.0 

0.994417 

0.994426 

2300.0 

0.994391 

0.994402 

2400.0 

0.994344 

0.994367 

2600.0 

0.994278 

0.994293 

2600.0 

0.994193 

Mmam 

2700.0 

0.994092 

0.994110 

TABLE  4.8  COMPARISON  OF  AVAILABILITY  FOR  VALUES  OF 
a  >  8.0,  |i  -  1460,  Rp  -  1  HOUR,  Rp  -  8  HOURS. 


Replaoement 
Age  (hours) 

Approadmate 

AvallabUity 

Exact 

Availability 

1600.0 

0.993961 

0.993963 

1600.0 

0.994189 

0.994192 

1700.0 

0.994373 

0.994376 

1800.0 

0.994617 

0.994521 

1900.0 

0.994628 

0.894632 

2000.0 

0.994709 

0.984714 

2100.0 

0.994764 

0.994770 

2200.0 

0.994796 

0.994802 

2300.0 

0.994806 

0.994812 

2400.0 

0.994796 

0.994803 

2600.0 

0.994767 

0.994777 

2600.0 

0.994722 

0.994733 

27000 

0.994662 

0.994674 

2800.0 

0.994687 

0.994601 

2900.0 

0.994498 

0.994514 

V.  POUCY  III  (MINIMAli  REPAIR  /  FAILURE  REPLACEMENT  / 
PREVENTIVE  REPLACEMENT) 


POLICY  II  assumes  that  each  time  a  failure  occurs  it  can  be  repaired  and 
the  system  is  restored  to  an  operational  state  without  changing  the  failure  rate 
of  the  system  as  a  whole.  Now  we  shall  be  more  realistic  and  consider  the 
possibility  of  a  me^or  failure  before  the  planned  age  replacement  time  t^*,  and 
that  failure  is  rectified  by  a  replacement  (Here  again  a  replacement  may  be 
an  overhaul  which  is  assumed  to  return  ihe  system  to  a  state  as  good  as  new). 
So  when  the  oystem  is  running,  two  types  of  failures  are  possible: 


•  Type  I  Mure,  denoted  by  1C.  This  Mure  is  corrected  by  minimal  repair; 
if  &e  Type  I  Mure  occurs  at  age  t,  the  failure  rate  just  after  correction 
is  h(t).  This  type  of  failure  is  usually  repaired  at  sea  by  the  ship's  staff 
using  the  support  elements  onboard. 

*  Type  II  Mure  denoted  by  2C.  This  failure  is  remedied  by  effective  system 
replacement;  the  Type  II  failure  is  followed  by  overhaul,  after  which  the 
failure  rate  is  h(0).  This  type  of  failure  is  beyond  the  ship's  staff 
capability  either  due  to  lack  of  expertise  or  unavailability  of  the  required 
support  elements,  and  the  ship  has  to  return  to  base  to  effect  repair 
ei^er  by  msgor  overhaul  or  replacement. 


The  failures  are  detected  immediately. 

It  is  assumed  that  the  costs  and  downtime  for  simple  preventive 
maintenance  actions  are  negligible. 


74 


First  let  us  assume  the  time  to  some  type  of  failure,  T,  has  an  IFR 
distribution  F(t)  with  a  failure  rate  denoted  by  h(t);  0  ^  t  is  the  system 
operating  time  since  last  overhaul.  Let  Y  be  the  time  of  a  type  2C  failure 
having  the  distribution  function  G(t). 

Let  Pi  be  the  probability  of  Type  1C  failure 
and  pg  be  the  probability  of  Type  2C  failure 
(pi  and  pg  could  dependent  on  time) :  Pi  Pa  <■  1«0 
Now  let  N(t)  denote  the  number  of  Type  1C  failures  that  occurs  in  time 
t,  where  t  is  measured  from  a  moment  of  replacement  or  overhaul.  By  the 
assumption,  the  expected  number  of  Type  1C  failures  in  time  t,  m(t)  ■  E[N(t)], 
satisfies  a  simple  differential  equation  obtained  as  follows: 
a(t4>dt)  a(t)(l  -  pih(t)dt)  -i-  {m(t)n>pxh(t)dt  •¥  o(dt)  (5.1) 
where  o(dt)  represents  higher  order  terms  which  are  ■  0 
rearranging  the  terms  we  get 

da(t) 

-  -  Pih(t) 

dt  (6.2) 


and  therefore 

«(0  -  Pi  J*  Kii) 


(5.3) 
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provided  minimal  repair  correction  times  are  assumed  to  be  negligible  enough 
BO  that  the  expected  downtime  from  that  source  is  dose  to  Ric'Zn(t),  where  Rjc 
is  the  mean  time  of  a  minimal  repair. 

Let  T  be  the  time  to  a  failure  (either  Type  1C  or  2C)  with  the  failure  rate 
h(t) 

Then  the  failure  rate  of  Type  I  failure  ■  Pih(t)  at  age  t 
and  the  flsilure  rate  of  Type  II  failure  ■  P2h(t)  at  age  t. 

Since 


and 


[1  -  m  ]  -  « 


(5.4) 


(1  -G(f)  1  - 

• 


HfX 


-A'*»** 


ti  -  0(0  1  -  (1  -  TO  (M) 

From  equation  (6.6),  it  can  be  seen  that  if  p, »  1,  then  the  distribution  of 

any  type  of  failure  is  the  same  as  the  distribution  of  a  Type  2C  failure  and  this 

implies  that  each  failure  is  a  msyor  failure  and  it  requires  a  replacement  / 

overhaul  action  and  this  is  the  case  as  in  Policy  I  (Age  Replacement)  as 

mentioned  in  Chapter  III. 

For  P2  ■  0,  this  impli^  i  that  the  probability  of  a  msgor  failure,  i.e  Type  2C 
is  0  and  each  failure  is  a  minimal  failure  and  the  failure  is  removed  by 
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minimal  repair  and  this  is  the  case  as  in  Policy  II  (Minimal  Repair  and 
Preventive  Replacement)  as  mentioned  in  Chapter  IV. 

For  the  case  when  0  <  pj  <  1,  then  a  failure  could  be  either 

•  Type  1C 

•  Type  2C 

This  is  Policy  III,  which  is  the  general  case  and  will  be  discussed  in  this 
chapter. 

This  idea  of  two  types  of  failure  has  been  studied  by  Beichelt  and  Fisher 
[Ref.  18].  They  derived  the  Reliability  functions  for  calculating  the  expected 
long  run  cost  rate  for  a  generalised  age*replacement  policy.  They  assumed  that 
maintenance  actions  take  only  negligible  times  which  in  reality  is  not  true, 
especially  when  a  msyor  failure  occurs  at  sea.  Besides  that,  in  the  military 
environment  we  are  often  interested  in  the  availability  of  the  system  as  our 
measure  of  effectiveness  so  in  this  poli<^  we  will  expand  the  cost  rate  model, 
but  also,  and  more  importantly,  formulate  a  model  to  maximize  the  availability 
of  the  system  subjected  to  two  types  of  failures. 

A.  COST  MODEL 

The  age  tp*  at  which  the  operating  system  is  replaced  or  overhauled 
depends  on  the  following  factors: 
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•  Failure  distribution  (IFR). 

*  Costs  of  itiinimal  repair,  preventive  replacement  and  failure  replacement. 

*  Downtime  of  preventve  replacement  and  failure  replacement. 

•  Measure  of  effectiveness: 

.  Minimize  costs. 

'  Maximize  availability. 

It  is  assumed  that  the  downtime  of  a  minimal  repair  is  negligible. 

Let  C^c  bo  ^0  cost  of  a  minimal  repair. 

Gao  be  the  cost  of  a  failure  replacement. 

Cap  be  the  cost  of  a  preventive  replacement. 

Rio  be  the  mean  time  of  a  minimal  repair  (T^po  IC  failure). 

Rac  be  the  mean  time  of  a  failure  replacement  (Type  2C  failure). 
Rap  be  the  mean  time  of  a  preventive  replacement,  tp.  be  the  planned 
replacement  age. 

N(tp)  be  the  niuuber  of  minimal  repairs  in  interval  (O.tp). 

y  be  the  observed  time  from  system  replacement  until  the  next 
Type  2C  failure. 

N(y)  be  the  number  of  minimal  repairs  in  the  period  (0,y),  where  y  is 
the  time  until  Type  2C  failure. 

Rae  >  Rap  >  R^o  and  Cpo  ^  Cpp  >  Cjc 

A  replacement  takes  place  either  at  time  tp  or  when  there  is  a  Type  2C 
failure. 
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1.  Case  1 


Y  >  t, 

I - X - X - X - 

0 

X— -minimal  rapair 
2.  Caae  2 


Y  <  t. 


X— —Animal  rapair 


A  cycle  is  completed  each  time  a  replacement  takes  place  and  the 
costs  C  incurred  in  a  cycle  is  given  by  the  total  costs  of  minimal  repairs  and 
the  cost  of  either  a  preventive  replacement  or  a  failure  replacement,  This  is 
given  by  the  follo^ng: 


C 


{  Cie  [  H(t,)  I  Y  >  t,  1  +  C,c 
{ 

<  [  H(y)  I  Y  51  t,  ]  +  Ca, 


w.p.  1  -  0(y) 

»r.p.  0(y)  (5.0) 


In  order  to  find  the  total  costs  of  minimal  repairs,  we  ought  to  know 
the  expected  number  of  minimal  repairs  and  this  is  given  below  for  both  the 
two  cases; 
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For  case  1,  Y  >  tp 


£[  I  F  >  g  -  /J'P,  h(t)  dt  (5.7) 

For  case  2,  Y  2  tp 
Let  N(y)  ■  n  ,  then 

«  ATCy)  I  F  i  -  X;  «  W(y)-n  |  Fs^ 

•-0  (5.8) 

•  n  f[A^(y)«ii ,  F 

“ft  J»[F  i  g 

thus 

£I  N(»)  I  r  <  1  -  «*)  *  1  *(0  *  (5.9) 

Then 

£[q  .  KCic  W(g  I  F  >  rg)  +  ][1-C?(gi  ,, ... 

KC^c  E{m  I  F  <  <g)  +  Cjc  l[0(g] 

The  length  of  the  cycle,  L,  depends  on  the  time  of  a  Type  2C  failure 
and  is  given  by: 


L  «{ 


Rgp  +  tp 

Rac  +  F 


y  >  tp 
Y  S  tp 


(5.11) 


Then  the  expected  length  of  the  cycle,  E[L],  is 

E[L]  -  +  ^o(g]  +  Afcg  (5.12) 

where  M(tp)  is  the  mean  age  of  the  system  at  replacement  /  overhaul  and  is 
given  by 
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Af(fp  -  t  git)  dt  +  tjil  -  0(g  ]  (5.13) 

which  can  be  simplified  to 

Mit;)  -  Q  [1  -  0(f)  ]  A  (5.14) 

From  equation  (3.1), 

*(C1 

C(tp)  - - 

l[X.] 

Then  the  optimal  value  of  tp*  is  the  value  that  minimizes  the  cost  rate 
fiinction  G(tp)  and  this  can  be  found  by  graphical  or  numerical  analysis.  [Ref. 
16]  highlights  that  an  finite  optimal  interval  tp*  exists  if  h(t)  is  monotone 
increasing  ftmction  and  C  >  0  where  C  is  given  by 

Cie 

C  ■  - 

(®ic  “  ®«f  “  Cie) 

8.  WeibuU  Example 

For  a  time  to  failure.  T,  following  a  Weibull  distribution  with  shape 
parameter  a  and  scale  parameter  A.  the  expected  total  cost  in  a  cycle.  E[C],  is 
given  as  follows: 
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(S.15) 


£[C]  .  ,-«*V  jc,cP,  (it/  t  C„  -  C„)  *  C„ 

,u-i  ^-ncur  j, 

The  integral  *'<>  has  no  dosed  form 


solution  and  the  IMSL  subroutine  DQDAG  was  used  to  solve  it  numerically. 
The  length  of  the  cyde  is  given  by 

E[L1  ■  dt  (5-l«) 

4.  Numerical  lUuatration 

Now  let  us  consider  an  example  with  these  data:  assuming  that  the 
downtime  for  a  minimal  repair  is  negligible: 

Cgp  m$  25,000 
*  $  37|600 

Cio  ■$  1,000 

Rap  ■  8  hours 
Rac  ■  16  hours 
Rjo  ■  1  hour 
a  ■  3.0 

X  ■  1/1390  hours 
Pi  ■  0.6 
Pa  ■  0.4 
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Figure  6.1  shows  the  plot  of  the  expected  cost  rate  function.  It  is 
observed  that  the  curve  is  fairly  flat  near  the  optimal  point  thus  giving 
flexibility  to  the  decision  maker.  The  optimal  replacement  age  is  tp  s  1888.64 
hours  at  a  cost  of  $22.03  per  hour. 

In  a  real  situation  P2  is  a  function  of  age  and  it  usually  increases 
with  age  i.e.  as  the  system  ages  the  probability  of  a  msgor  failure  approaches 
1.0.  For  our  case  we  have  taken  the  probability  to  be  constant  and  this  again 


is  to  simplify  our  computation. 


Figure  6.1  Plot  of  the  Coat  Rate  Function. 
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TABLE  5.1  VALUES  OF  COST/HOUR  FOR  VARIOUS 
REPLACEMENT  AGE. 


Replacement  age 
tp*  (hours) 


1100.0 


1200.0 


1300.0 


1400.0 


1600.0 


1600.0 


1700.0 


1800.0 


1900.0 


2100.0 


2200.0 


2300.0 


2400.0 


2600.0 


Coet 

($  /  hour) 


27.67 


26.01 


24.76 


23.83 


23.16 


22.67 


22.36 


22.16 


22.06 


22.03 


22.06 


22.13 


22.22 


22.32 


22.42 


22.62 
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B.  AVAILABILITY  MODEL 


4 


i 


Now  let  us  take  tlie  (Point)  Availability  as  our  measure  of  effectiveness. 
When  the  ship  is  at  sea  there  are  ixsually  two  types  of  failures;  failixres  which 
can  be  rectified  onboard,  each  of  which  individually  does  not  effect  the 
availability  of  the  system  drastically,  and  other  failures  which  are  not 
repairable  onboard  so  the  ship  has  to  return  to  port  to  effect  repairs;  this 
causes  the  availability  to  be  more  severely  degraded.  First  we  shall  assume 
that  the  downtime  for  a  minimal  repair  is  negligible.  Then  using  the  equation 
(3.8),  the  availability  function  is  given  by 


(5.17) 


where  M(tp)  is  the  mean  operational  time  in  a  cycle  and  it  is  given  by 

M(ip  ■  [1  -  (7(1)  ]  A  (5.18) 

and  this  is  similar  to  the  availability  function  in  Policy  I  in  Chapter  III. 

Again  the  availability  function  in  equation  (6.17)  is  an  approximation 
because  we  do  not  take  into  account  the  downtime  during  a  minimal  repair, 
which  in  reality  is  not  faithful  to  reality,  but  should  be  reasonably  accurate. 
If  we  do  take  the  downtime  of  the  minimal  repair  into  account  we  can 


obtain  an  exact  solution  as  follows: 


86 


Using  the  same  principle  as  in  equation  (4.15), 

Let 

Y  ■  Time  of  a  Mayor  Failure 

{ 1  if  the  system  is  up  at  time  t 

X(t) «  { 

(  0  if  the  system  is  dovm  at  time  t 

P  «  P{Pailtire  is  Type  1C) 

ai(t)  -  P{Y  >  t,  X(t)  -  1) 

ao(t)  -  P{Y  >  t,  X(t)  -  0) 

ai(t)  +  ao(t)-P{Y>t}  -l-G(t) 

{X(t)  ifY>t 
Z(t) «  { 

{0  ifY.<t 

h(t)  ■  the  failure  rate  of  Type  1C  failure,  which  is  pih(t) 
p(t)  H  the  repair  rate  of  Type  1C  failure,  and  for  the  Exponential 
distribution,  p(t)  ■  Rio'^ 

Now, 

ai(t-t-dt)  «  ax(t){l  -  h(t)<it}  *  ao(t)H(t)dt  *  o(dt)  (6.19) 

ao(b+dt)  -  ao(t)(l  -  ^(t)dt)  -I-  ai(t)h(t)dt  -i-  o(dt)  (6.20) 

Initial  conditions;  ax(0)  ■  1  ;  ae(0)  ■  0 
Equation  (6.19)  becomes 

d(ax) 

- -  -  h(t)ax(t)  +  ao(t)n(t)  (6.21) 

dt 
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-  ^l(t)«o(t)  +  ai(t)li(t)p  (5-22) 

dt 

Since  ao(t)  ■  {1  -  G(t))  -  ai(t),  Equation  (6.21)  becomes 

- -  -  b(t)«i(t)  +  -  0(t)]  -  ai<t)>  (6.28) 

dt 

which  can  be  simplified  to 

- -  -  {h{t)  +  ^l(t)^ax{t)  +  p(t){l  -  G(t)}  (6.24) 

dt 

The  equation, 

-  +  {h(t)  +  ^t(t))  ai(t)  «  0  /  with  ai(0)  -  1 

dt 

has  the  solution 

ai(t)  ■  BXPC-  l(t)} 

where 

fi(f)  ■  * 

Therefore  the  equation  (6.24)  has  the  solution 

fliCO  -  (1  +  {1  -  0(5)}  ds]  (6.26) 

and  equation  (6.26)  can  now  be  written  as 
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(6.26) 


P,(0  .  €  <*! 

and  for  the  Weibidl  distribution  with  shape  parameter  a  and  scale  parameter 


X,  and  substituting  the  values  for  h(t)  and  p(t)  in  equation  (6.26),  the 
probability  that  the  system  is  up  at  time  t  is  given  by 


Pi(t)  ”  e 


But  the  integral 


-  Pi(u^  *  - -p-: 


has  no  closed-form 


expression  in  terms  of  elementary  tabulated  functions.  Therefore  we  have  to 
evaluate  it  numerically, 

Since  the  length  of  the  cycle,  L  is 


{  Raf  +  tp 
•( 

(  Rao  + 


Y  >  tp 

Y  i  tp 


then  using  the  same  principle  as  that  used  in  [Hef.  16]  the  average  availability 
over  a  cycle  of  length  L  and  hence  in  the  long  nm  is  given  by 


il(rp) 


(5.28) 


Eiq)ecmd  length  of  a  Cycle 
But  the  Expected  length  of  a  cycle  is  given  by  equation  (6.12).  Therefore 
the  Availability,  A(tp),  is  given  by 
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(6.29) 


ilOp 


_ /;*  O.M  <fc _ 

J^ll  -  0(«pi  ♦  JljclO(»pi  +  /o'  1  -  <* 


To  study  the  above  expression  the  IMSL  subroutine  DQDAG  was  used  to 
compute  the  Availability  function  a^Ct)  for  various  values  of  t.  Then  we  used 
the  IMSL  subroutine  DGSINT  [Ref.  13]  to  compute  cubic  spline  interpolant  to 
the  set  of  data  points  obtained  earlier  (values  of  t  and  the  corresponding  values 
of  a^Ct) ).  Finally  we  used  the  IMSL  subroutine  DCSITG  [Ref.  13]  to  evaluate 
the  integral  of  the  cubic  spline  for  various  values  of  the  replacement  age  tp. 
From  this  we  can  compute  the  optimal  replacement  age  tp*  and  the  maximum 
availability.  The  values  obtained  by  using  this  'exact*  solution  was  compared 
with  the  approximation  as  in  equation  (6.17).  These  values  are  tabulated 
below.  It  is  observed  that  both  the  approximate  and  exact  solutions  are  the 
same  to  three  decimal  places,  and  the  approximation  gives  an  accurate 
solution,  However  this  is  true  for  small  values  of  Rjo  and  Rgc*  As  the  values  of 
Rio  and  R^c  are  increased  the  accuracy  of  the  present  approximation  method 
degrades. 

Table  6.2  shows  the  values  of  availability  for  a  b  3.0,  p  ■>  1390,  Ric  ■  1 
hr,  Rgo  ■  16  hrs,  Rgp  ■  8  hrs.  The  optimal  replacement  age  for  the  two  methods 
and  the  maximum  availability  are: 


89 


•  approximate  solution  1628  hours  0.991715 

•  exact  solution  1496  hours  0.991567 

If  we  use  the  approximate  solution  as  our  reference  for  the  replacement  action, 
the  exact  availability  is  0.991666,  which  is  extremely  dose  to  the  exact 
solution,  and  much  more  easily  obtained. 

Table  6.3  shows  the  values  of  availability  for  a  «  3.0,  p  ■  139C.  Rjc  *  8 
hr,  Rao  B  24  hrs,  Rgp  ■  8  hrs.  The  optimal  replacement  age  for  the  two  methods 
and  the  maximum  availability  are: 

•  approximate  solution  1201  hours  0.989796 

•  exact  solution  1113  hours  0.989070 

If  we  use  the  approximate  solution  as  our  reference  for  the  replacement  action, 
the  exact  availability  is  0.9890 16  which  is  extremely  close  to  the  exact  solution, 
Table  6.4  shows  the  values  of  availability  for  a  «  3.0,  p  s  1390,  Rjc  ^  8 
hr,  Rac  s  72  hrs,  Rjp  *  8  hrs.  The  optimal  replacement  age  for  the  two  methods 
and  the  maximiun  availability  are: 

•  approximate  solution  762  hours  0.984143 

•  exact  solution  736  hoturs  0.983836 

If  we  use  the  approximate  solution  as  our  reference  for  the  replacement  action, 
the  exact  availability  is  0.983827  which  is  very  close  to  the  exact  solution. 
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Table  6.6  shows  the  values  of  availability  for  a  s  3.5,  \x  s  1390,  Ric  »  8 
hr,  Rjic  =  24  hrs,  R^p  =  8  hrs.  The  optimal  replacement  age  for  the  two  methods 
and  the  maximum  availability  are: 

•  approximate  solution  1147  hours  0.990119 

•  exact  solution  1076  hours  0.989619 

If  we  use  the  approximate  solution  as  our  reference  for  the  replacement  action, 
the  exact  availability  is  0.989471  which  is  again  very  dose  to  the  exact 
solution. 

Table  6.6  shows  the  values  of  availability  for  a  »  2.6,  p  «  1390,  R}c  ■  8 
hr,  Rjo  ■  24  hrs,  Rjp  m  8  hrs.  The  optimal  replacement  age  for  the  two  methods 
and  the  maximum  availability  are: 

•  approximate  solution  1317  hours  0.989493 

•  exact  solution  1200  houia  0.988594 

If  we  use  the  approximate  solution  as  our  reference  for  the  replacement  action, 
the  exact  availability  is  0.988536  which  is  once  again  very  dose  to  the  exact 
'  solution. 

»  From  the  above  results  we  can  condude  that  the  approximate  solution  is 

a  good  approximation  for  planning  the  replacement  or  overhaul  actions  of  a 
system  in  order  to  maximize  the  availability  of  the  system.  From  oiu*  earlier 
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discussion  we  found  that  the  availability  is  sensitive  to  the  shape  parameter 
a;  Figure  6.2  shows  the  availability  (exact  solution)  for  various  values  of  alpha. 

If  the  shape  parameter  a  was  in  fact  3.6  and  in  our  estimation  we  used 
a  B  2.6,  then  we  would  have  lost  an  availability  of  0.000626,  which  is 
extremely  small.  This  shows  that  the  parameters  need  not  be  estimated  very 
accurately  to  achieve  good  results. 

Figure  6.2  shows  that  for  low  values  of  alpha  such  as  2.6,  the  availability 
function  is  rather  flat  and  there  is  more  flexibility  in  determining  the 
replacement  age  that  is,  the  replacement  interval  at  1200  hours  or  1500  hours 
gives  about  the  same  availability  on  the  long  run.  However  this  is  not  true  for 
higher  values  of  alpha  such  as  3.6  where  the  availability  fhnction  falls  quite 
rapidly  on  both  sides  of  the  optimal  replacement  age  tp*.  At  the  availability  of 
about  0.9883,  the  replacement  age  is  1475  hours  and  it  is  insensitive  to  the 
value  of  alpha. 

Now  we  shall  also  look  at  the  effects  of  the  scale  parameter  k  on  the 
availability.  Table  6.7  shows  the  values  of  availability  for  a  s  3.0,  p  ■  1360,  R^c 
»  8  hr,  Rgc  ■  24  hrs,  Rap  «  8  hrs.  The  optimal  replacement  age  for  the  two 
methods  and  the  maximiim  availability  are: 

•  approximate  solution  1167  hours  0.989496 

•  exact  solution  1081  hours  0.988760 
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Figure  5.2  The  Plot  of  TSxaot’  Availability  for  Various  Values  of  the 

Shape  Parameter  Alpha. 


If  we  use  the  approximate  solution  as  our  reference  for  the  replacement 
action,  the  exact  availability  is  0.988693  which  is  very  dose  to  the  exact 
solution. 

Now  if  the  actual  scale  parameter  was  1390  hours  and  it  was  estimated 
to  be  1360  hours,  then  we  would  have  lost  an  availability  of  0.00032  which  can 
be  considered  very  small. 

Table  6.8  shows  the  values  of  availability  for  a  s  3.0,  p  b  1460,  Ric  «  8 
hr,  Rac  »  24  hrs,  Rgp  *  8  hrs.  The  optimal  replacement  age  for  the  two  methods 
and  the  maximiim  availability  are: 
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TABLE  5.2  COMPARISON  OF  AVAILABIUTY  FOR  VALUES  OF 
01  ■  8«0y  |i  ■  1890|  Rjo  ■  1  HR)  Rfc  *  HRS)  R^p  *  8  HRS 


Replacement 
Age  (houn) 

Approximate 

Availability 

Exact 

Availability 

1320.0 

0.991692 

0.991474 

1344.0 

0.991621 

0.991499 

1368.0 

0.991646 

0.991620 

1392.0 

0.991666 

0.991637 

1416.0 

0.991683 

0.991649 

1440.0 

0.991696 

0.991669 

1464.0 

0.991706 

0.991664 

1488.0 

0.991711 

0.991667 

1612.0 

0.991716 

0.991666 

1636.0 

0.991716 

0.991663 

1660.0 

0.991713 

0.991667 

1684.0 

0.991709 

0.991649 

1608.0 

0.991702 

0.991638 

1632.0 

0.991693 

0.991626 

1666.0 

0.991682 

0.991611 

TABLE  5.8  COMPARISON  OF  AVAILABILITY  FOR  VALUES  OF 
Q  B  8.0(  )i  ■  1890(  Ric  B  8  HRSf  Ric  **  HRS^  Rgp  ■  8  HRS 


Replacement 

Age  Gioim) 

Approximate 

Availability 

Exact 

Availability 

1018.0 

0.989670 

0.988997 

1042.0 

0.989629 

0.989030 

1066.0 

0.989678 

0.989063 

1090.0 

0.989718 

0.989066 

1114.0 

0.989749 

0.989070 

1138.0 

0.989772 

0.989066 

1162.0 

0.989787 

0.989062 

1186.0 

0.989794 

0.989032 

1210.0 

0.989796 

0.989004 

1234.0 

0.989789 

0.988970 

1268.0 

0.989778 

0.988929 

1282.0 

0.989760 

0.988882 

1306.0 

0.989737 

0.988830 

1330.0 

0.989709 

0.988773 

1364.0 

0.989677 

0.988710 
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TABLE  6.4  COMPARISON  OF  AVAILABILITY  FOR  VALUES  OF 
a  «  3.0,  |i  ■  1890,  Ric  -  8  HRS,  P  c  -  72  HRS,  R«p  -  8  HRS 


Replacement 
Age  (houra) 

Approximate 

Availability 

Exact 

Availability 

520.0 

0.982290 

0.982140 

644.0 

0.982706 

0.982642 

668.0 

0.983067 

0.982877 

692.0 

0.983349 

0.983164 

616.0 

0.983688 

0.983376 

640.0 

0.983779 

0.983650 

664.0 

0.983926 

0.983680 

688.0 

0.984033 

0.983768 

712.0 

0.984102 

0.983819 

736.0 

0.984138 

0.983836 

760.0 

0.984141 

0.983818 

784.0 

0.984116 

0.983772 

808.0 

0.984061 

0.983697 

832.0 

0.983981 

0.983596 

866.0 

0.983877 

0.983470 
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TABLE  S.5  COMPARISON  OF  AVAILABILITY  FOR  VALUES  OF 
a  -  8.6,  VI  -  1890,  R^o  -  8  HRS,  R^  -  24  HRS,  Rv  -  8  HRS 


Replacement 
Age  (houra) 

Approximate 

AvailabiUty 

9440 

0.989738 

968.0 

0.989828 

992.0 

0.989905 

1016.0 

0.989969 

1040.0 

0.990021 

1064.0 

0.990061 

1088.0 

0.990091 

1112.0 

0.990110 

1136.0 

0.990119 

1160.0 

0.990119 

1184.0 

0.990110 

1208.0 

0.990092 

1232.0 

0.990067 

1266.0 

0.990034 

1280.0 

0.989994 

Exact 

Availability 

0.989330 

0.989394 

0.989446 

0.989482 

0.989506 

0.989617 

0.989617 

0.989606 

0.989484 

0.989462 

0.989411 

0.989361 

0.989302 

0.989235 

0.989160 
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TABLE  6.6  COMPARISON  OF  AVAILABH^ITY  FOR  VALUES  OF 
(X  ■  2.6)  |i  ■  1800)  ■  8  HRS)  B^c  *  HRS)  R^p  ■  8  HRS 


Replacement 
Age  (houre) 

Approsdmate 

Availability 

Exact 

Availability 

1120.0 

0.989330 

0.988661 

1144.0 

0.989370 

0.98S679 

1168.0 

0.989404 

0.988689 

1192.0 

0.989431 

0.988694 

1216.0 

0.989454 

0.988693 

1240.0 

0.989471 

0.988687 

1264.0 

0.989482 

0.988576 

1288.0 

0.989490 

0.988660 

1312.0 

0.989493 

0.988640 

1336.0 

0.989492 

0.988616 

1360.0 

0.989487 

0.988487 

1384.0 

0.989478 

0.988466 

1408.0 

0.989466 

0.988421 

1432.0 

0.989461 

0.988383 

1466.0 

0.989433 

0.988342 

TABLE  6.7  COMPARISON  OF  AVAILABILITY  FOR  VALUES  OF 
(X  ■  8.0)  pi  w  1860y  Rjo  ■  8  BRS)  R^p  ■  24  HRSi  R^p  ^  8  HRS 


Replacement 
Age  (houra) 

Approximate 

Availability 

Exact 

Availability 

948.0 

0.989132 

0.988588 

972.0 

0.989214 

0.988643 

996.0 

0.989284 

0.988687 

1020.0 

0.989343 

0.988718 

1044.0 

0.989392 

0.988738 

1068.0 

0.989430 

0.98S748 

1092.0 

0.989469 

0.988749 

1116.0 

0.989480 

0.988740 

1140.0 

0.989492 

0.988722 

1164.0 

0.989496 

0.988697 

1188.0 

0.989494 

0.988664 

1212.0 

0.989484 

0.988623 

1236.0 

0.989468 

0.988677 

1260.0 

0.989446 

0.988623 

1284.0 

Mmam 

0.988464 

TABLE  6.8  COMPARISON  OF  AVAILABILITY  FOR  VALUES  OF 
a  -  8.0,  |i  -  1450,  Ric  -  8  HRS,  Rk  -  24  HRS,  R^  -  8  HRS. 


Replacement 
Age  (hours) 

Approximate 

Availabillly 

Exact 

Availability 

0.989812 

0.989322 

1024.0 

0.989892 

0.989379 

1048.0 

0.989961 

0.989426 

1072.0 

0.990021 

0.989461 

1096.0 

0.990072 

0.989488 

1120.0 

0.990114 

0.989606 

1144.0 

0.990148 

0.989616 

1168.0 

0.990175 

0.989516 

1192.0 

0.990194 

0.989611 

1216.0 

0.990207 

0.989497 

1240.0 

0.990213 

0.989478 

1264.0 

0.990213 

0.989462 

1288.0 

0.990208 

0.989420 

1312.0 

0.990197 

mmM 

1336.0 

0.990181 

0.989340 

•  approximate  solution  1263  hours  0.990214 

•  exact  solution  1161  hours  0.989616 

If  we  use  the  approximate  solution  as  our  reference  for  the  replacement 
action,  the  exact  availability  is  0.989464  which  is  again  very  close  to  the  exact 
solution. 

Now  if  the  actual  scale  parameter  was  1450  hours  and  it  was  estimated 
to  be  1360  hours,  then  we  would  have  lost  an  availability  of  0.000446  which 
once  again  can  be  considered  very  small.  This  shows  that  the  scale  parameter 
does  not  effect  the  availability  drastically  and  a  close  estimate  is  sufficient. 

C.  MAXIMIZEAVAILABIUTYSUBJECTTOBUDGETCONSTRAINT 
The  availability  and  cost  are  very  important  measures  of  effectiveness, 
vt/e  would  like  to  have  as  many  resources  as  possible  to  maximize  the 
availability  of  a  system,  however  in  reality  we  are  often  limited  by  budget 
constraints.  So  we  would  like  to  achieve  cost  effectiveness,  that  is  we  would 
like  to 

maximize  Effectiveness  Level 
subject  to  Budget  ^  B 

The  effectiveness  level  which  is  commonly  used  is  the  availability.  So  in 
our  case  we  would  like  to 
Maximize  A(tp) 

Subject  to  C(tp)  ^  B. 
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We  have  solved  this  problem  graphically.  Figure  6.3  shows  the  plots  of  the 
availability  and  the  cost  rate  functions  us  formulated  in  Section  B.  It  is 
observed  that  the  availability  is  maximum  at  0.991567  at  an  optimal 
replacement  interval  tp*  of  1496  hours  while  the  cost  is  minimum  at 
$22.03/hour  at  an  optimal  replacement  interval  tp*  of  1889  hours. 

Now  if  we  have  a  budget  of  not  exceeding  $22.25/hour,  then  the  maximum 
availability  that  can  be  obtained  is  C  9915  with  an  optimal  replacement 


interval  tp*  of  1650  hours. 


Function. 

The  other  measture  of  effectiveness  that  is  also  used  is  the  mission 
reliability  and  this  again  will  give  anoUier  optimal  replacement  age,  so  now  we 
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have  a  multiple  independent  conilicting  criterion  and  it  is  up  to  the  Decision 
Maker  to  decide  which  of  the  measitre  of  effectiveness  is  vital  and  fits  the 
scenario  very  well. 
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VL  CONCLUSION  AND  RECOMMENDATIONS 


The  fighting  efiTectiveneas  and  operational  readiness  of  a  ship  depends 
largely  on  the  operational  availability  of  her  equipment  and  systems.  If  we  are 
not  constrained  by  the  budget  then  we  can  expend  the  resources  necessary  to 
achieve  the  desired  availability,  such  as  incorporating  more  redundancies  or 
carrying  out  "premature"  replacements  end  overhauls.  However  in  reality  this 
is  not  the  case,  and  we  are  always  limited  by  the  available  budget,  and  we 
would  like  to  maximize  system  availability  subject  to  budget  constraints. 

At  sea,  equipment  and  systems  are  exposed  to  various  environment  and 
unfavourable  conditions.  As  such,  they  are  subject  to  stochastic  failure  and 
deterioration.  However,  with  timely  maintenance  actions  as  discussed  in  the 
various  policies  in  the  thesis,  we  can  minimize  catastrophic  and  unexpected 
failures,  enabling  us  to  achieve  the  desired  measures  of  effectiveness. 

Therefore,  based  on  the  optimal  maintenance  polides,  we  can  carry  out 
replacement  actions  or  complete  overhaul  of  the  equipments  and  systems  at 
the  base  during  the  stand-off  periods  so  that  when  the  ship  is  at  sea  we  can 
minimize  loss  of  availability  due  to  failures  and  maximize  our  successes  in  the  , 

operational  missions. 

Based  on  the  policies,  we  can  also  carry  adequate  spares.  This  is 
particularly  important  for  long  missions  or  when  the  need  to  be  prepositioned 
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in  forward  operating  areas.  On  the  other  hand,  although  we  can  carry  spares 
onboard,  sometimes  other  support  elements  are  not  available  to  rectify  the 
defects. 

In  our  analysis,  we  have  assumed  that  we  know  the  failure  distribution. 
For  our  case  we  concentrated  on  the  WeibuU  distribution.  The  two  important 
parameters  which  are  usually  estimated  are  the  shape  parameter,  alpha,  and 
the  scale  parameter,  lambda.  From  the  analysis  it  is  observed  that  these 
parameters  need  not  be  estimated  very  accurately.  A  slight  variation  in  the 
values  of  the  parameters  do  not  affect  the  long  run  availability  and  cost  rate 
functions  drastically.  Of  the  two,  the  shape  parameter  needs  to  be  estimated 
more  accurately. 

In  the  formulation  of  the  availability  functions,  we  simplified  the 
computation  by  assuming  the  downtime  for  a  minimal  repair  to  be  negligible. 
When  compared  with  the  exact  solution,  taking  into  accoimt  the  downtime  for 
a  minimal  repair,  the  results  obtained  by  the  approximation  method  gave 
extremely  accurate  results.  Many  of  the  fimctions  and  integrals  that  were 
formulated  did  not  have  closed  form  expressions  in  terms  of  elementary 
tabulated  functions.  However  IMSL  subroutines  were  available  in  the  Math 
library  at  the  main  frame  at  the  Naval  Postgraduate  School  and  these 
subroutines  expeditiously  computed  the  integrals  very  accurately. 

The  expected  downtimes,  as  taken  in  our  analysis,  are  practical  figures 
assuming  all  the  support  elements  are  readily  avalilable  when  reqxiired.  That 
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is  the  reason  for  the  high  values  of  the  availability.  In  reality  we  often  have  to 
wait  for  spare  parts  and  support  elements,  which  sometimes  have  long  lead 
times. 

In  the  thesis  we  also  studied  the  effects  of  simple  preventive  maintenance 
actions  on  the  measiires  of  effectiveness.  It  is  observed  that  simple  preventive 
maintenance  actions  do  not  restore  the  system  to  a  condition  "as  good  as  new" 
but  the  maintenance  actions  can  enhance  or  improve  the  reliability  of  the 
system  by  a  certain  factor  which  decreases  as  the  system  ages.  However,  a 
methodology  for  chacterizing  the  effective  age  reduction  remains  to  be 
developed.  In  our  anal}rsis  we  have  assumed  that  the  system  improves  by  a 
certain  factor  on  completion  of  each  preventive  maintenance  action. 

Most  planned  maintenance  systems  usually  adopt  maintenance  efforts 
based  on  calendar  time  (weeks,  months)  or  running  hours  of  systems  or 
equipment,  but  from  the  analysis  it  is  observed  that  for  systems  that  have  a 
"wear-out"  life  distribution  we  shall  have  to  successively  resort  to  decreasing 
maintenance  intervals  if  we  are  going  to  maintain  the  systems  above  some 
minimum  reliability  level. 

This  subject  can  be  expanded  farther  by  future  research.  The  following 
areas  are  recommended: 

•  Carry  out  similar  analysis,  especially  for  the  availability  function  for  the 
three  policies  when  the  underlying  life  distribution  F  comes  from  a  family 
of  Qomma  Distribution. 
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*  Incorporate  imperfect  repair  into  the  model. 

*  A  methodology  for  chacterizing  the  improvement  of  the  reliability  of  the 
system  on  completion  of  preventive  maintenance  needs  to  be  investigated. 

*  A  system  may  consists  of  various  sub-systems.  The  maintenance  actions 
for  the  sub-systems  has  to  be  coordinated  so  that  instead  of  just  taking 
into  account  the  availability  of  the  sub-systems  individually,  the  whole 
system  has  to  be  considered  bearing  in  mind  of  some  dependence  on  the 
supporting  elements  associated  with  the  sub-systems. 


* 
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APPENDIX  A.  COMPUTER  PROGRAM  AND  OUTPUT 


PROGRAM  PM 

C  THIS  PROGRAM  IS  FORMULATED  TO  COMPUTE  THE  AVAILABILITY, 

C  RELIABILITY  AND  COST  FUNCTIONS  FOR  VARIOUS  POLICIES  OP 
C  INTEREST.  THE  PROGRAM  DETERMINES  THE  OPTIMAL  REPLACEMENT 
C  OR  OVERHAUL  INTERVAIB  FOR  THE  SYSTEM.  THE  OPTIMAL 
C  REPLACEMENT  INTERVAL  DEPENDS  ON  THE  MEASURES  OF 
C  EFFECTIVENESS  DESIRED  AND  IHE  RELEVANT  TYPE  OF  POLICY 
C  APPLICABIJB  TO  THE  SYSTEM  UNDER  STUDY.  THE  PROGRAM 
C  UTILIZES  THE  IMSL  SUBROUTINES  FROM  THE  MATH  UBRARY 
C  AVAILABLE  AT  THE  MAIN  FRAME  AT  THE  NAVAL  POST  GRADUATE 
C  SCHOOL.  IT  IS  ASSUMED  THAT  THE  TIME  TO  FAILURE  FOLLOWS 
C  THE  WEIBULL  DISTRIBUTION  WITH  SHAPE  PARAMETER  ALPHA  AND 
C  SCALE  PARAMETER  LAMBDA. 


PARAMETER  STATEMENT 
INTEGER  NDATA,NTINT 
PARAMETER  (NDATA  -  81,  NTINT  -  20) 


VARIABLES  DECLARATION 

REALMS  A,  B,  F,  RESULT,  ERRABS,  ERRREL,  ERREST,  ALPHA, 
&  LAMBDA,  MU,  C2C,  CIC,  C2P,  R2P,  R2C,  RIP,  RIC, 

&  Al,  AA,  XDATA(NDATA),  YDATA(NDATA), 

St  BREAK(NDATA),  CSCOEF(4,NDATA),  AVAIL,  AVAILl, 

&  PI,  P2,  30LN,  SF.  CDF,  LENGTH,  G,  H,  ANSWER, 

&  DCSITG,  COST,  COSMIN,  ANS,  D,  REL,  E,  EE,  FF, 

&  TD2P,  TDTl,  TDT2,  SUMDIP,  SUMD2C,  STIME,  TL, 

&  CDPTL,  TXNTINT),  REL,  R,  TDT,  UPTIME,  TDIP, 

&  SUM,  SUMl,  X,  XX.  LL,  UL,  STEP 

INTEGER  IRULE,POL,P,SPMREQ,NN 

C  VARIABLES  DEFINITION 
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OOOOOOOOOOOOOOOOOOOOOOOODOOOOOOOOOOOOOOOOOOOOC3 


A . LOWER  LIMIT  OF  INTEGRATION  USED  IN 

ARGUMENT  OF  THE  IMSL  SUBROUTINE  DQDAG  AND 
IMSL  FUNCTION  DCSITG 

B . UPPER  LIMIT  OF  INTEGRATION  USED  IN 

ARGUMENT  OF  THE  IMSL  SUBROUTINE  DQDAG  AND 
IMSL  FUNCTION  DCSITG 
E^E,P,PP,G,H..PUNCTIONS  TO  BE  INTEGRATED 

RESULT . ESTIMATE  OF  THE  INTEGRAL  FROM  A  TO  B  OF 

THE  FUNCTIONS  E,BE,F,PF 

ERRABS . ABSOLUTE  ACCURACY  DESIRED  AS  THE  INPUT 

ARGUMENT  OP  THE  IMSL  SUBROUTINES  DQDAG 

ERRREL . RELATIVE  ACCURACY  DESIRED  AS  THE  INPUT 

ARGUMENT  OF  THE  IMSL  SUBROUTINES  DQDAG 

mULE . CHOICE  OF  QUADRATURE  RULE  (SEE  APPENDDC 

2) 

ALPHA . SHAPE  PARAMETER  OP  WEIBULL  DISTRIBUTION 

LAMBDA . SCALE  PARAMETER  OF  WEIBULL  DISTRIBUTION 

MU . RECIPROCAL  OF  LAMBDA 

C2C . EXPECTED  COST  OF  A  FAILURE  REPLACEMENT 

CIC . EXPECTED  COST  OF  A  MINIMAL  REPAIR 

C2P . EXPECTED  COST  OP  A  PREVENTIVE  REPLACEMENT 

R2C . EXPECTED  DOWNTIME  OP  A  FAILURE 

REPLACEMENT 

RIC . EXPECTED  DOWNTIME  OF  A  MINIMAL  REPAIR 

R2P . EXPECTED  DOWNTIME  OP  A  PREVENTIVE 

REPLACEMENT 

A1 . COMPUTE  EXPONENTIAL  EXPRESSION  (POLICY 

II) 

AA . AVAILABIUTY  AT  TIME  T  (POLICY  II) 

NDATA . NUMBER  OF  DATA  POINTS  FOR  COMPUTING  CUBIC 

SPLINE  INTERPOLANT.  INPUT  ARGUMENT  FOR 
IMSL  SUBROUTINE  DCSINT 

XDATA(  ) . ARRAY  OP  LENGTH  NDATA  CONTAININO  THE  DATA 

POINTS  ABSCISSAS.  INPUT  ARGUMENT  FOR  IMSL 
SUBROUTINE  DCSINT 

YDATA(  ) . ARRAY  OP  LENGTH  NDATA  CONTAININO  THE  DATA 

POINTS  ORDINATES.  INPUT  ARGUMENT  FOR  IMSL 
SUBROUTINE  DCSINT 

BREAK(  ) . ARRAY  OF  LENGTH  NDATA  CONTAINING  THE 

BREAKPOINTS  FOR  THE  PIECEWISE  CUBIC 
REPRESENTATION.  OUTPUT  ARGUMENT  OP  THE 
IMSL  SUBROUTINE  DCSINT 

CSCOEP . MATRIX  OP  SIZE  4  BY  NDATA  CONTAININO  THE 

LOCAL  COEFFICIENTS  OF  THE  CUBIC  PIECES. 

OUTPUT  ARGUMENT  OP  THE  IMSL  SUBROUTINE 
DCSINT 


C  AVAIL . APPROXIMATE  AVAILABILITY 

C  ANSWER . EXACT  AVAILABILITY 

C  PI . PROBABILITY  OF  TYPE  I  (MINIMAL)  FAILURE 

C  P2 . PROBABILITY  OF  TYPE  II  (MAJOR)  FAILURE 

C  SF . SURVIVAL  FUNCTION 

C  CDF . DISTRIBUTION  FUNCTION 

C  LENGTH . LENGTH  OF  A  CYCLE 

C  COST,  COSMIN...COST  RATE  FUNCTION 

C  D . MISSION  DURATION 

C  REL . RELIABILITY 

C  TD2P,TDT2 . TOTAL  DOWNTIME  OP  PREVENTIVE  REPLACEMENT 

C  TD1P,TDT1 . TOTAL  DOWNTIME  OP  PREVENTIVE  MAINTENANCE 

C  SUMDIP . TOTAL  DOWNTIME  WHEN  TIME  OP  FAILURE  IS 

C  AFTER  PREVENTIVE  REPLACEMENT 

C  8UMD2C . TOTAL  DOWNTIME  WHEN  TIME  OF  FAILURE  IS 

C  BEFORE  PREVENTIVE  REPLACEMENT 

C  STIME . TOTAL  DOWNTIME  IN  A  CYCLE 

C  UPTIME . TOTAL  UPTIME  IN  A  CYCLE 

C  LL . LOWER  LIMIT  OF  REPLACEMENT  INTERVAL 

C  UL . UPPER  UMIT  OP  REPLACEMENT  INTERVAL 

C  STEP . STEP  SIZE  OR  INCREMENT 

C  POL . POUCY  YO  BE  EVALUATED 

C  SPMREQ . SIMPLE  PREVENTIVE  MAINTENANCE  REQUIREMENT 

C  EXTERNAL  SUBROUTINES  AND  FUNCTIONS 
EXTERNAL  E.EE,F,FF,G,HJ)QDAO,DCSINT,DCSITG 


C  COMMON  BLOCKS 

COMMON /FT/ R1C,B 
COMMON  /GT/  P2 
COMMON  /IT/  PI 
COMMON  /HT/  LAMBDAALPHA 


C  INITIALIZATION 

A  -  O.ODO 
ERRABS  -  O.ODO 
ERRREL  -  O.OOOOIDO 
mULE  -2 
PREL  -l.ODO 
SUMDIP  -  O.ODO 
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SUMD2C  =  O.ODO 
STIME  e  O.ODO 
SUM  « O.ODO 


C  Prompt  the  user  for  the  policy  to  be  evaluated 
POUCY I - AGE  REPLACEMENT 

POLICY  n - AGE  REPLACEMENT  WITH  MINIMAL  REPAIR 

POLICY  m  —  AGE  REPLACEMENT  WITH  TWO  TYPES  OF 
FAILURES.  'PYPE  I  MINIMAL  FAILURE,  TYPE 
II  MAJOR  FAILURE 

PRINT*, TINTER  THE  NUMBER  OP  THE  POLICY  TO  BE  EVALUATED 
READ*J»OL 

PRINT*,TOLICY  -  ’,POL 
GOTO  (100,200,300),POL 

00  CONTINUE 

PRINT*,’DO  YOU  WANl'  TO  INCLUDE  SIMPLE  PREVENTIVE 
&  MAINTENANCE-ENTER  1  FOR  NO  AND  2  FOR  YES’ 
READ«,SPMREQ 
GOTO  (126,1S0).SPMREQ 


POLICY  I  (AGE  REPLACEMENT-NO  SIMPLE  PREVENTIVE 
MAINTENANCE) 


126  CONTINUE 

PRINT*,ENTER  THE  ESTIMATED  VALUES  FOR  THE  SHAPE 
&  PARAMETER  ALPHA  AND  THE  SCALE  PARAMETER  MU  ’ 
READ*,ALPHA,MU 
PRINT*, ’ALPHA  «i’,ALPHA,’MU  «  ’,MU 

PRINT*, ’ENTER  ESTIMATED  VALUES  FOR  THE  EXPECTED  COST  OP 
&  PREVENTIVE  REPLACEMENT  FOLLOWED  BY  THE  EXPECTED 
&  COST  OF  FAILURE  REPLACEMENT 

READ*,C2P,C2C 

PRINT*, TEXPECTED  COST  OP  A  PREVENTIVE  REPLACEMENT  - 
&  ’,C2P 

&  , ’EXPECTED  COST  OP  A  FAILURE  REPLACEMENT  « 

&  ’,C2C 

PRINT*, ’ENTER  VALUES  FOR  THE  EXPECTED  DOWNTIME  OF 


111 


&  PREVENTIVE  REPLACEMENT  FOI.LOWED  BY  EXPECTED 

&  DOWNTIME  OF  A  FAILURE  REPLACEMENT 

READ'*32P,R2C 

print*;expected  downtime  of  a  preventive  replacement 

&  ’412P 

&  ,  'EXPECTED  downtime  OF  A  FAILURE  REPLACEMENT  = 

&  *32C 

PRINT*, 'ENTER  THE  VALUE  OP  MISSION  DURATION  TIME  IN 
&  HOURS’ 

READ*4) 

PRINT*,'  MISSION  DURATION  TIME  *  D,  'HOURS’ 

PRINT*,’ENTER  THE  VALUE  FOR  LOWER  LIMIT,  UPPER  UMIT, 

&  STEP  SIZE  ’ 

READ*,LL,UL,STEP 

PRINT*,’L0WER  limit  .  ’,LL,’UPPER  LIMIT  -  ’,UL,'STBP  - 
&  ’.STEP 

WRITE(S0,101)ALPHA,MU,C2P,C2G,R2P,R2C,D 

101  P0RMAT(16X,T0UCY’,3X,T,3X,’(AGE  REPLACEMENT)' 

&  /16X,30(’-’)y/, 

&  10X,’ALPHA  «  ’,P3.1,17X,  'MU  r  ',P6.1/, 

&  10X,’C2P  -  $’,P7.1,12X,  ’C2C  -  $’,P7.i;, 

&10X,’R2P  -’,P4.1,’H0URS’,11X,’R2C  -  ’,P4.1,' 

&  HOURS 

&  V.IOX.'MISSION  DURATION  ■  ’,P4.1,'  HOURSV/, 

&  6X,’REPLACEMENT.6X,’COST  RATE’,6X,’AVAILABILITr,6X, 

&  'MISSION  RELIABILITY’,/, 

&  6X,’AGE  (HOURS)’,6X,’($/HOUR)’,10X,’(A)’,  17X,’(R)'/, 

&  5X,11(’-’),  6X,8(’.’),  0X,12(’-’),  6X.19(’-')) 

LAMBDA  «  l.ODO/MU 
DO  11 1  -  LL,UL,STEP 
B-DBLE(I) 

CALL  DQDAG  (E,A,B,ERRABS,ERRREL,IRULE,RESULT,ERREST) 
SF«EXP(-((LAMBDA*B)**ALPHA)) 

CDF-l.ODO-SP 

COST-((C2P*SF)+(C2C*CDP)y((B2P*SP)+(R2C*CDP)+RESULT) 

AVAIL-RESULT/((R2P*SP>f(R2C*CDF)+RESULT) 

REL»  DEXP(-((LAMBDA*(B+D))**ALPHA))  / 

&  DEXP(-<(LAMBDA*B)**ALPHA)) 

WRITE(30,102)B,COST>^VAIL,REL 

102  PORMAT(6X,P7. 1,9X,P8.3,8X,P8.6,12X,P8.6) 
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o  ooo 


11  CONTINUE 


GOTO  99 


POUCY  I  (AGE  REPLACEMENT  -  WITH  SIMPLE  PREVENTIVE 
MAINTENANCE) 


160  CONTINUE 

PRINT*, ENTER  THE  ESTIMATED  VALUES  FOR  THE  SHAPE 

&  PARAMETER  ALPHA  AND  THE  SCALE  PARAMETER  MU ' 

READ*,ALPHA,MU 

PRINT*, 'ALPHA  -’,ALPHA,’MU  « ’,MU 

PRINT*,ENTER  VALUES  FOR  THE  EXPECTED  DOWNTIME  OF 
&  PREVENTIVE  REPLACEMENf  FOLLOWED  BY  EXPECTED 

&  DOWNTIME  OF  A  FAILURE  REPLACEMENT 
READ*J12P,R2C 

PRINT*, 'EXPECTED  DOWNTIME  OP  A  PREVENTIVE  REPLACEMENT  - 
&  ',R2P 

&  ,  'EXPECTED  DOWNTIME  OF  A  FAILURE  REPLACEMENT  - 

&  '32C 

PRINT*, 'ENTER  THE  EXPECTED  DOWNTIME  FOR  SIMPLE 
&  PREVENTIVE  MAINTENANCE  AND  THE  IMPROVEMENT 
&  FACTOR  ON  COMPLETION  OP  THE  MAINTENANCE' 
READ*31P»R 

PRINT*,THE  DOWNTIME  FOR  SIMPLE  PREVENTIVE  MAINTENANCE 
&  -  '31P 

&  ,’THE  IMPROVEMENT  FACTOR  *  ',R 

PRINT*,'ENTER  THE  VALUE  FOR  LOWER  LIMIT,  UPPER  LIMIT, 

&  STEP  SIZE' 

READ*,LL,UL,STEP 

PRINT*,'LOWER  UMIT  -  ’,LL,'UPPER  UMIT  -  ',UL,'STEP  - 
&  ',STEP 

WRITE(30,103)ALPHA,MU4WP32C,R1P,R 
103  FORMA'r(6X,TOLICY’,2X,'r42X,’(AQE  REPLACEMENT  WITH 
&  SIMPLE',1X,  'PREVENTIVE  MAINTENANCE)',/,8X, 

&  63  ('-')  // 

&  10X,'ALPHA  «  ',F3.1,17X,'MU  «  ',F6.1,/, 

&  10X,'R2P  =  ',P4.1,'  HOURS', 11X,'R2C  «  ',F4.1,' 
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&  HOURS’, y.lOX, 111?  =  ’,F4.1,’  HOURM2X,’R  = 

&  ’,F4.1,//, 

&  6X,’NO.  OP  INTERVAL’, SX.’REPIACEMENT  AGE  (HOURS)’, 6X, 
&’AVAILABIUTY’y,6X,16(’-’),6X,2S(’-’),6X,12(’-’)) 

LAMBDA  »  l.ODO/MU 

T(l)  -  (1.0D0/LAMBDA)*((-(LOG(0.76D0)))** 

&  (l.ODQ/ALPHA)) 

CDF  -  l.ODO  -  EXP(-((LAMBDA*T(1))**ALPHA)) 

SP  « l.ODO-CDP 
B  «T(1) 

DO  10  I «  1,NTINT 
IP  (I  .GT.  1)  THEN 
T(I)-T(I.1)+(((1.0D0-R)**(M))*T(1)) 

END  IF 

10  CONTINUE 

DO  20  N  -  LL,UL,STEP 
IP  (N  .EQ.  1)  THEN 
TDT  -  (R2P*SP)+(R2C*CDF) 

A  -  O.ODO 

CALL  DQDAQ(FF,A,B,ERRABS,ERRREL,IRULE, RESULT, 

&  ERREST) 

UPTIME  -  (ALPHA*(LAMBDA**ALPHA)*RESULT)+(SP*T(1)) 
EI^E 

DO  30  J  «  1,N 
IF  (J  .EQ.  1)  THEN 
TL  -  O.ODO 
ELSE 

TL  -  T(1)^T(J).T(J-1)) 

END  IF 

CDPTL  ■  l-EXI^-CdiAMBDA-^D^^ALPHA)) 

REL  -  1-(CDF.CDFTL) 

PREL  .  PREL*REL 
30  CONTINUE 

TDIP  -  (N-1)*R1P*PREL 
TD2P  -  R2P*PREL 
TDTl  ■  TD1P+TD2P 
PREL  «  l.ODO 

DO  40  K  «  2;S 
DO  60  L  -  1,K-1 
IP  (L  .EQ.  1)  THEN 
TL  «  0.0D0 
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ELSE 

TL  s  T(lHT(L)-T(Lrl)) 

END  IP 

CDPTL  =  l.EXP(.((LAMBDA*TL)**ALPHA)) 

REL  «  l.(CDF-CDPTL) 

PREL  -  PREL*REL 
60  CONTINUE 

TL  -  T(1MT(K)-T(K-1)) 

CDPTL  *  l.EXP(.((LAMBDA*TL)**ALPHA)) 

FT  «  CDF-CDFTl. 

SUMDIP  »  SUMDIP+^K-D^RIP^PREL-^PT) 

SUMD2C  =  SUMD2C+(R2C*PREL*PT) 

PREL  -  l.ODO 
40  CONTINUE 

TDT2  «  SUMD1P+SUMD2C+(R2C*CDF) 

TDT-  TDT1+TDT2 
PREL  «  l.ODO 
DO  60  M  -2.N 
DO  70  P  -  l,M.l 
IP  (P  .EQ.  1)  THEN 
TL  -  O.ODO 
ELSE 

TL  -  T(l).(T(P).T(P.l)) 

END  IF 

CDPTL  -  1-EXP(-((LAMBDA*TL)**ALPHA)) 

REL  -  l.(CDP-CDPTL) 

PREL  -  PREL'^REL 
70  CONTINUE 

A  «  T(1MT(M).T(M-1)) 

CALL  DQDAQ{FF.A,B.ERRABS,ERRREL,IRULE, RESULT, 
&  ERREST) 

ANSWER  «  ALPHA*(LAMBDA**ALPHA)*RESULT 
TL  -  T(l)-(T(M).T(M.l)) 

CDPTL  -  l.EXP(-<(LAMBDA*TL)**ALPHA)) 

REL  -  l.(CDF.CDFTL) 

ANSWER  -  ANSWER+(REL*(T(M).T(M-1))) 

STIME  -  STIME  +(PREL*ANSWER) 

PREL  -  l.ODO 
60  CONTINUE 


A  -  O.ODO 

CALL  DQDAG  (FFUV,B,FHRABS,ERRREL,IRULE, 
&  RESULT, ERREST) 
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ANSWER  »  (ALPHA*(LAMBDA**ALPHA)*RESULT)  + 

$  (SP*T(1)) 

UPTIME  »  STIME-i-  ANSWER 
END  IF 

AVAIL  -  UPTIME/(UPTIME+TDT) 

WRITE(30,104IN,T(N)AVAIL 
20  CONTINUE 

104  FORMAT(10X.I3,19X.F6.1,17X,F8.6} 

GOTO  99 

C  POUCY  n  (AGE  REPLACEMENT  -  WITH  MINIMAL  REPAIR) 

200  CONTINUE 

PRINT^.’ENTER  THE  ESTIMATED  VALUES  FOR  THE  SHAPE 
&  PARAMETER  ALPHA  AND  THE  SCALE  PARAMETER  MU  ’ 

R£AD*U^lU.MU 
PRINT*, ’ALPHA  «’,ALPHA,’MU  =  ’,MU 

PRINT*, ’ENTER  ESTIMATED  VAI^UES  FOR  THE  EXPECTED  COST  OP 
&  PREVENTIVE  REPLACEMENT  FOLLOWED  BY  THE  EXPECTED 
&  COST  OF  MINIMAL  REPAIR’ 

READ*,C2P,C1C 

PRINT*, TEXPECTED  COST  OP  A  PREVENTIVE  REPLACEMENT  - 
&  ’,C2P 

&  .’EXPECTED  COST  OF  A  MINIMAL  REPAIR 
&  ’,C1C 

PRINT*, ENTER  VALUES  FOR  THE  EXPECTED  DO^'NTIME  OF 
&  PREVENTIVE  REPLACEMENT  FOLLOWED  BY  EXPECTED 
&  DOWNTIME  OF  A  MINIMAL  REPAIR’ 

READ*,R2P,R1C 

PRINT*, EXPECTED  DOWNTIME  OP  A  PREVENTIVE  REPLACEMENT  » 
&  ’,R2P 

&  ,  EXPECTED  DOWNTIME  OP  A  MINIMAL  REPAIR 

&  ’.RIC 


PRINT*, ENTER  THE  VALUE  FOR  LOWER  LIMIT,  UPPER  LIMIT, 
&  STEP  SIZE’ 

READ*,LL,UL,STEP 

PRINT*, ’LOWER  LIMIT  =  ’.LL.’UPPER  LIMIT  «  ’,UL,’STEP  - 
&  ’3TEP 


oooo  ^  ooo 


LAMBDA  -  l.ODO/MU 

PRINT*, DO  YOU  WANT  TO  INCLUDE  SIMPLE  PREVENTIVE 
&  MAINTENANCE-TYPE  1  FOR  NO  AND  2  FOR  YES’ 
READ*,SPMREQ 

GOTO  (225,260)3PMREQ 


POUCY II  (AGE  REPLACEMENT  -  WITH  MINIMAL  REPAIR  WITHOUT 
SIMPLE  PREVENTIVE  MAINTENANCE) 


25  CONTINUE 

WRITE(80,105)ALPHA,MU,C2P.C1C32P31C 
06  FORMAT(9X,*POUCY’,2X.’ir,2X,’(AGE  REPLACEMENT  WITH 
&  MINIMAL*, 1X,’REPAIR)’,/,9X50(’-V4 
&  10X,’ALPHA  -  ’,P8.1,17X,*MU  -  ’,P6.1,/, 

&  10X,’C2P  »  $*F7.1,12X,*C1C  «  $’,P7.1,/, 

&  10X,’R2P  -  ’,P4.1,’  H0URS’,11X,’R1C  «  ’,P4.1,* 

&  HOURS*;/, 6X,’REPLACEMENT,0X.’COSr,lOX,’APPROXIMATE*, 

&  8X,’EXACr;,  5X,’AGE  HOURS)*, 5X,*$/HOUR*,9X,  * 

&  AVAILABILITY*,  4X,  *AVAILABILITY’;,6X,11(*-’), 

&  5X,7(’-*),  fiX.12(*.*),4X,12(*.*)) 

DO  12  I »  0,4000,60 
B-DBLE(I) 

CALL  DQDAG  (EEAB,ERRABS,ERRREL,IRULE  JiESULT,ERREST) 
A1  «DEXP(-(((LAMBDA*B)**ALPHA)4<B/R1C))) 

AA  -  A1+(RESULT*(1.0DO/R1C) ) 

WRITE(18,*)B,AA 
L2  CONTINUE 

REWIND  (18) 

DO  22  J  »  1,  NDATA 
READ(18,*)XDATA(J),YDATA(J) 

22  CONTINUE 

CALL  DCSINT(NDATA,XDATA,YDATA,BRBAK,CSCOEF) 


Calculate  the  integral  of  the  spline 
approximation. 
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OOQO 


DO  32  K  «  LL,UL,STEP 
B  .  DBLE(K) 

NINTV  -  NDATA  - 1 

ANSWER  -  DCSITG(A,B,NINTV3REAK.CSC0EF) 
ANSWER  .  ANSWER/(B-i-R2P) 

AVAIL  -  (B^R1C*((LAMBDA*B)**ALPHA)))/(B+R2P) 
AVAILl »  B/(B-i-R2P-»-(RlC*((LAMBDA*B)**ALPHA))) 
GOST  « ((ClC*((LAMBDA*B)*«ALPHA)HC2PV(B-t-R2P) 
WRITE  (30406)  B,C0ST^VAILANSW£R 
32  CONTINUE 

106  FORMAT(7X>F6.1,6X,F9.440X,F8.6, 8X38.6) 

GOTO  99 


I 


POUCY  n  (AGE  REPLACEMENT  -  WITH  MINIMAL  REPAIR  AND 
SIMPLE  PREVENTIVE  MAINTENANCE) 


260  CONTINUE 

PRINT*, TINTER  THE  EXPECTED  CCBT  OF  A  SIMPLE  PREVENTIVE 
&  MAINTENANCE  FOLLOWED  BY  THE  IMPROVEMENT  IN  THE 
&  AGE  IN  HOURS  ACHIEVED  ON  COMPLETION  OF  THE 
&  MAINTENANCE* 

READ*,C1P3X 

PRINT*, ’EXPECTED  COST  OP  SIMPLE  PREVENTIVE  MAINTENANCE 
&  -  ’,C1P, 

&  IMPROVEMENT  IN  AGE(HOURS)  -*XX 

PRINT*,’ENTER  THE  NUMBER  OP  SIMPLE  PREVENTIVE 
&  MAINTENANCE  ACTIONS  TO  BE  TAKEN  IN  A  CYCLE  ’ 
READ*3N 

PRINT*,’NUMBER  OF  SIMPLE  PREVENTIVE  MAINTENANCE  ACTIONS 
&  -’3N 

PRINT*, ’ENTER  THE  VALUE  FOR  LOWER  LIMIT,  UPPER  UMIT, 

&  STEP  SIZE’ 

READ*,LL,UL,STEP 

PRINT*,’LOWER  UMIT  -  ’,LL,’UPPBR  LIMIT  -  ’,UL,’STEP  - 
&  ’,STEP 

WRITE(30,107)ALPHA,MU,C2P,C1C,C1P3X,NN 
107  PORMAT(9X,’POLICY’,2X,’ir,2X,’(AGE  REPLACEMENT  WITH 
&  MINIMAL  REPAIR  AND  SIMPLE  PM  ) 

&  ’,/,9X,66(’-’)y/, 

&  10X,’ALPHA  -  ’,P3.1,17X,’MU  -  ’,F6.i;, 
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&  10X,’C2P  =  $’^7.1,12X,’C1C  =  $’,F7.V, 

&  10X,’C1P  «  $’^7.1,12X,’XX  =  ’^6.1/  HOURSV, 

&  lOX’NO.  OP  SIMPLE  PM  = 

&  6X,’REPLACEMENT',6X,’COSr;,  5X,’AGE  (HOURS)’, 6X,’ 

&  $/HOUR  ’  /6X,11(’.’),6X,7(’.’)) 

DO  6  J  -  LL,UL3TEP 
B  >  REAK J) 

x-xx 

IP  ((X-B)  .GB.  0)  THEN 
X.B 
END  IP 

DO  7 1  -  OW-1 

SUMl «  C1C*(LAMBDA**ALPHA)*(((B+((B.X)*I))**ALPHA) 
&  -  (((B-X)<'*ALPHA)*(I**ALPHA))) 

SUM  -  SUM -»■  SUMl 
7  CONTINUE 

COST  -  (SUM  +  ((NN-1)*C1P)  +C2P)  /  ((NN*B)  +R2P) 
LENGTH  «  B*REAL(NN) 

WRITE<30.108)LENGTH,OOST 
108  PORMAT(5X,P6.1,5X,F0.4) 

SUM  «  0.0 
6  CONTINUB 

GO  TO  99 


POLICY  HI  (AGE  REPLACEMENT  -  WITH  TWO  TYPES  OP  FAILURE) 


300  CONTINUE 

PRINT*, ’ENTER  THE  ESTIMATED  VALUES  POR  THE  SHAPE 
&  PARAMETER  ALPHA  AND  THE  SCALE  PARAMETER  MU  ’ 
READ*,ALPHA,MU 
PRINT*,’ALPHA  -’,ALPHA,’MU  -  ’,MU 

PRINT*,’ENTER  ESTIMATED  VALUES  B’OR  THE  EXPECTED  COST  OP 
St  PREVENTIVE  REPLACEMENT  B'OLLOWED  BY  THE  EXPECTED 
&  COST  OP  FAILURE  REPLACEMENT,  AND  COST  OP  MINIMAL 
&  REPAIR* 

READ‘>',C2P,C2C,C1C 

PRINT*,’EXPECTED  COST  OF  A  PREVENTIVE  REPLACEMENT  « 
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&  ’,C2P 

&  , 'EXPECTED  COST  OF  A  FAILURE  REPLACEMENT  - 

&  ’  C2C 

&  , 'EXPECTED  COST  OF  A  MINIMAL  REPAIR 

&  ’.CIC 

PRINT*-, 'ENTER  VALUES  FOR  THE  EXPECTED  DOWNTIME  OP 
&  PREVENTIVE  REPLACEMENT  FOLLOWED  BY  EXPECTED 

&  DOWNTIME  OF  A  FAILURE  REPLACEMENT  AND  DOWNTIME 

&  OF  MINIMAL  REPAIR’ 

READ-»32P,R2C.R1C 

PRINT*, 'EXPECTED  DOWNTIME  OF  A  PREVENTIVE  REPLACEMENT  - 
&  '32P 

&  ,  'EXPECTED  DOWNTIME  OF  A  FAILURE  REPLACEMENT  « 

&  ’32C 

&  ,  'EXPECTED  DOWNTIME  OF  A  MINIMAL  REPAIR 

&  ’4WC 

PRINT*, 'ENTER  THE  VALUE  FOR  LOWER  LIMIT,  UPPER  LIMIT, 

&  STEP  SIZE  ’ 

READ*.LL,UL,STEP 

PRINT*, ^OWER  LIMIT  .  ’,LL,'UPPER  LIMIT  .  ',UL,'STEP  - 
&  ',STBP 

PRINT*,’ENTER  VALUES  FOR  THE  PROBABIUTY  OF  MINIMAL 
&  FAILURE  AND  MAJOR  FAILURE' 

RBAD*J»1,P2 

PRn^*,THE  PROBABILITY  OF  MINIMAL  FAILURE  -  Tl, 

&  AND  PROBABILITY  OF  MAJOR  FAILURE  -  'P2 

PRINT*,’ENTER  THE  VALUE  FOR  LOWER  LIMIT,  UPPER  UMIT, 

&  STEP  SIZE  ’ 

READ*,LL,UL,STEP 

PRINT*,'LOWER  UMIT  «■  ’,LL,’UPPER  UMIT  «  ’,UL,’STEP  « 

&  ’,STEP 

LAMBDA  «  l.ODO/MU 

WRITE(30,109)AU>HA,MU,C2P,C2C,C1C,R2P,R2CJ11C,P1,P2 
109  PORMAT(9X,'POLICY,2X.’nr,2X,’(AGE  REPLACEMENT  WITH 
&  TWO  TYPES  OF  FAILURE’/, 9X,66(’.’V/, 

&  $  6X,’ALPHA  -  ',FS.1,11X,’MU  -  ’,F6.1/, 

&  6X,’C2P  -  $’,F7.1,6X,’C2C  *  $’,P7.1,7X.'C1C  - 
&  $’,F7.1/,  6X,’R2P  -’,F4.1,*HOURS’,6X,'R2C  - 
&  ’,F4.1,’ HOURS’,6X, 'RIC  -  ’,F4.1,’ HOURS’/,6X, 

&  TROBABILITY  OF  TYPE  I  FAILURE  -’,F3.1/,SX, 
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&  'PROBABILITY  OF  TYPE  II  FAILURE  = 

&  ,5X,'REPLACEMENr,6X,'COST',10X,'APPROXIMATE',8X, 

&  'EXACry,  5X,'A0E  (HOURS)’  ,5X,’$/HOUR’»9X,’ 

&  AVAILABILITY’,4X,’AVAILABILnY;,6X,ll('-’), 

&  4X,8(’-’),8X42(’-’)>  &  4X,12('-')) 

DO  13  1-0,4000,60 
B-DBLE(I) 

C  CALL  IMSL  SUBROUTINE  DQDAO 

CALL  DQDAG  (F,A3.ERRABS,ERRREL,IRULE, RESULT, ERREST) 

A1 »  DEXP(-<(((LAMBDA*B)**ALPHA)*Pl>f(B/RlC))) 

AA  -  A1-KRESULT*(1.0DQ/R1C) ) 

C  OUTPUT  THE  AVAILABILITY  AT  TIME  T  TO  PILE  13 

WRITE(1S,29)BM 
29  F0RMAT(1X,F6.1,16X,F8.6) 

18  CONTINUE 

REWIND  (13) 

DO  43  J  -  1,  NDATA 
READ(18,*)XDATA(J),YDATA(J) 

43  CONTINUE 

CALL  DCSINT(NDATA,XDATA.YDATA.BREAK,CSCOEF) 


Calculate  the  integral  of  the  ipline 
approximation. 

DO  33  K  -  LL,UL,STEP 
B  -  DBLE(K) 

NINTV  -  NDATA  - 1 

CALL  DQDAG  (G,A,B,ERRABS,ERRREL,IRULE,SOLN,ERREST) 
SP  -  DEXP(^P2*((LAMBDA*B)**ALPHA))) 

CDF  -  l.ODO-SF 

LENGTH  -  (R2P'*SP)+(R2C*CDF>+SOLN 
AVAILl  -  SOLN/LENGTH 

ANSWER  -  DCSITG(A,B,NINTV3REAK,CSC0EF) 

AVAIL  -  ANSWER/LENGTH 

CALL  DQDAG  (H,A,B,ERRABS,ERRREL,IRULE  ANS,ERREST) 
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COST  -  (C10*SP*(P1*(LAMBDA*B)**ALPHA))  +  (C2P*SP)  + 
&  (C2C*CDP)  +  (ALPHA*P1*P2*C1C'»(LAMBDA*'* 

&  (2.0D0*ALPHA))  *ANS) 

COSMIN  -  COST  /  LENGTH 

WRITE  (30,110)  B.COSMINU^VAILIU^VAIL 
38  CONTINUE 

110PORMAT(7X,F6.1,6X,F9.6.10X.P8.6.8X,F8.6) 

99  STOP 
END 


DOUBLE  PRECISION  FUNCTION  E(T) 
REAL'^8  LAMBDAUajPHA,T 
COMMON  /HT/  LAMBDA,ALPHA 
E-DEXP(-((LAMBDA  *  T)**ALPHA)) 
RETURN 
END 


DOUBLE  PRECISION  FUNCTION  BE(T) 

REAL«8  LAMBDA,ALPHA,TJ11C,B 
COMMON  /HT/  LAMBDA,ALPHA 
COMMON /FT/ R1C,B 

EE-DEXP(-(((LAMBDA  *  B)««ALPHA)-((LAMBDA«T)**ALPHA)  -i- 
&  (B/RIC) 

&-<T/RlC))) 

RETURN 

END 


DOUBLE  PRECISION  FUNCTION  P(T) 

REAL*8  LAMBDA,ALPHA,T31C,B.P1  J>2 
COMMON/PT/  R1C3 
COMMON /HT/  LAMBDA,ALPHA 
COMMON/IT/ PI 
COMMON /OT/  P2 

P-DEXP(-((((LAMBDA  *  B)**ALPHA)*P1)  - 
&  (((LAMBDA*T)**ALPHA)*P1)+(B/R1C)  -(T/RIC)  + 
&  (((LAMBDA*T)*'*ALPHA)*P2) )) 

RETURN 

END 


* 
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DOUBLE  PRECISION  FUNCTION  FF(T) 

REALMS  LAMBDA, ALPHA, T 

COMMON  /HT/  LAMBDA,ALPHA 

PF-(T**ALPHA)*(DEXP(-((LAMBDA  *  T)**ALPHA))) 

RETURN 

END 


DOUBLE  PRECISION  FUNCTION  G(T) 

REAL'^8  LAMBDA,ALPHA,T,P2 

COMMON/GT/P2 

COMMON /HT/  LAMBDA, ALPHA 

G-DEXP(.(P2*((LAMBDA*T)**ALPHA))) 

RETURN 

END 


DOUBLE  PRECISION  FUNCTION  H(T) 
REALMS  LAMBDA,ALPHA,T,P2 
COMMON/C3T/P2 
COMMON /HT/  LAMBDA,ALPHA 
H-(T**((2.0D0*ALPHA).1.0D0))  * 

&  DEXP(<P2*((LAMBDA*T)**ALPHA))) 

RETURN 

END 
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ALPHA  -  3.0 
C2P  -  $26000.0 
R2P  -  8.0  HOURS 
MISSION  DURATION 


MU  «  1390.0 
C2C  $37500.0 
R2C  16.0  HOURS 
24.0  HOURS 


REPLACEMENT  COST  RATE  AVAILABILITY  MISSION  RELIABILITY 

AflEiHQVRg)  msm  (A) _ (R) 


900.0 

32.781 

0.988396 

0.977947 

960.0 

31.900 

0.988565 

0.976493 

1000.0 

31.178 

0.988681 

0.972916 

1060.0 

30.692 

0.988762 

0.970214 

1100.0 

30.123 

0.988786 

0.967391 

1160.0 

29.764 

0.988786 

0.964447 

1200.0 

29.472 

0.988761 

0.961383 

1260.0 

29.263 

0.988718 

0.968200 

1300.0 

29.115 

0.988649 

0.964899 

1860.0 

29.020 

0.988671 

0.951482 

1400.0 

28.968 

0.988483 

0.947961 

1460.0 

28.961 

0.988389 

0.944806 

1600.0 

28.962 

0.988290 

0.940649 

1660.0 

28.996 

0.988191 

0.936681 

1600.0 

29.043 

0.988003 

0.932704 
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1  (AGE  REPLACEMENT  WITH  SIMPLE  PREVENTIVE  MAINTENANCE 


ALPHA 

-  3.0 

MU  «  1890.0 

R2P  » 

8.0  HOURS 

R2C  B  48.0  HOURS 

RIP  « 

1.0  HOUR 

o 

n 

NO.  OP  INTERVAL 

REPLACEMENT  AfiE..(H0UR8)  AYAILAPiLITy 

1 

817.6 

0.979421 

2 

1743.4 

0.982197 

8 

2486.7 

0.983306 

4 

3156.6 

0.983691 

6 

8767.7 

0.983686 

6 

4299.6 

0.983480 

7 

4787.1 

0.988344 

8 

6226.0 

0.983206 

9 

5621.0 

0.983076 

10 

5976.6 

0.982866 

11 

6296.5 

0.982848 

12 

6684.4 

0.982761 

13 

6848.6 

0.982664 

U 

7076.8 

0.982687 

16 

7286.7 

0.982618 
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POLICY  II  (AGE  REPLACEMENT  WITH  MINIMAL  REPAIRl 


ALPHA 

»  3.0 

MU  « 

C2P  « 

$25000.0 

CIC  « 

R2P  - 

8.0  HOURS 

RIC  s 

REPLACEMENT 

COST 

APPROXIMATE 

AOEIHgURS) 

msm 

AVAILABIMTY 

1200.0 

40.6693 

0.992845 

1300.0 

42.6667 

0.993268 

1400.0 

44.9682 

0.993693 

1500.0 

47.8289 

0.993862 

1600.0 

51.1154 

0.9940/6 

1700.0 

54.8016 

0.994246 

1800.0 

58.8682 

0.994374 

1900.0 

63.2987 

0.994469 

2000.0 

68.0807 

0.994632 

2100.0 

73.2039 

0.994569 

2200.0 

76.6598 

0.994661 

2300.0 

84.4416 

0.994671 

2400.0 

90.6432 

0.994640 

2500.0 

96.9601 

0.994490 

2600.0 

103.6880 

0.994423 

126 


1390.0 
$37600.0 
1.0  HOUHS 


EXACT 

avami^ty 

0.992847 

0.993261 

0.993595 

0.993865 

0.994060 

0.994249 

0.994379 

0.994474 

0.994539 

0.994576 

0.994589 

0.994580 

0.994551 

0.994503 

0.994437 


POLICY  III  (AGE  REPLACEMENT  WITH  TWO  TYPES  OF  FAILURE 


ALPHA  -  3.0  MU  »  1890.0 

C2P  »  $26000.0  C2C  »  $37600.0  CIC  -  $  1000.0 

R2P  «  8.0  HOURS  R2C  -  24.0  HOURS  RIC  »  8.0  HOURS 

PROBABILITY  OP  TYPE  I  FAILURE  «  0.6 
PROBABIUTY  OF  TYPE  II  FAILURE  »  0.4 


REPLACEMENT 

COST 

APPROXIMATE 

EXACT 

AQE  (HOURS) 

AVAILABILITY 

AVAILABILITY 

1018.0 

27.30690 

0.989670 

0.988997 

1042.0 

26.88668 

0.989629 

0.989080 

1066.0 

26.49167 

0.989678 

0.989068 

1090.0 

26.12318 

0.989718 

0.989066 

1114.0 

26.77884 

0.989749 

0.989070 

1138.0 

26.46711 

0.989772 

0.989066 

1162.0 

26.16666 

0.989787 

0.989062 

1186.0 

24.87626 

0.989794 

0.989082 

1210.0 

24.61474 

0.989796 

0.989004 

1234.0 

24.37107 

0.989789 

0.988970 

1268.0 

24.14424 

0.989778 

0.988929 

1282.0 

23.93333 

0.989760 

0.988882 

1306.0 

23.78747 

0.989737 

0.988830 

1330.0 

23.66683 

0.989709 

0.988773 

1364.0 

23.38767 

0.989677 

0.988710 

1378.0 

23.23226 

0.989639 

0.988643 
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APPENDIX  B.  DERIVATION  OF  COST  RATE 
AND  AVAILABILITy  FUNCTIONS 


Let  T  be  the  time  to  failure 


<  C,  T  >  t, 

C  ■  ( 

{  C,  Tit, 

Then  I  [C]  «  C,  P[T>t,]  *  Cf 

-  Cj  P[1  -  P<tp)l  +  C,  PlP(t,)l 

(  +  t,  T  >  t, 

L  -  ( 

{  T  T 

Thenl[L]  -  I  [I.  I  T  >  t,]  P[P  >  t,l  -f  1[L  I  V  £  t,]  P[T  i  t^l 


Now, 

PtT  i  t 


Then, 


T  £  ty]  ■  P[T  1  t/  T  S  t^] 

PI*  S  t,l 

(  P(t) 

(  - 

-  (  Pit,) 

<  1 


it  t  i  t, 
if  t  >  tp 


£[r  I  7  <  g 


Jo  /ny 


Therefore 

1[L]  -  (IV  +  VII  -  f(Vl  +  *»  IP  (VI  +  /o'  A 
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■  [1  -  F(tp)]  -f  »,  I»{tp)]  +  t,  [1  -  F(tp)J  + 

■  »,  [i  -  r(t,)i  +  I,  +  t,  -  t,  F(t,))  + 

«A0  * 

Th«  integral  t  f(]t)  dt  can  be  simplified  by  using 

J  u  dv  ■  uv  “  J  V  du 
then  f  XO  A  -  t,  F  (tp)  -  m  dt 

now  l[L]  -  Up  [1  -  F(tp)]  +  »,  [F(tp)l  +  tp  -  m  dt 

-  Up  II  -  F(tp)l  +  »,  IF(tp)l  +  f^^dt  - 

■  Up  [1  -  F<tp)l  +  *1  E»(t,)l  +  /J'  [1  “  F(fi]  dt 

-  Up  [1  -  F(tp)l  +  F,  IF(tp)l  +  ll(tp> 

where  M(tp)  ■  J*'  [1  -  i»l(0]  dt 

■  IC] 

Now  C(tp)  ■  - 

I[L] 

_ Cp  F[1  -  F(tp)l  »  C,  P[F(tp)l _ 

Rp  [1  -  F(tp)l  +  »,  IF(tp)l  +  JJ'  [1  -  F(0]  dt 
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Setting  dC(tp) 

- ■  0 

dtp 

we  get  C*  £(V  [1  -  Fm  [1  " 

C.  t»(t,)l  -  C,  [F(V1  +  C,  IF(t,)3»  -  C,  IF{t,)l* 

-  Cp  -  cv  (F(V1  +  *<V  [»A  -  »pC,l 
rearranging  the  terms  we  get 

Iq  [i  -  ^01  ^  <C*  -  C,)  -  F(tp)  {1  -  F(V )  (C*  -  CV) 

C,  (1  -  Ftt,)}  +  £(t,)  [»*C,  -  B^C,] 

£(V 

since  h  (tp)  ■ - 

1  -  F(t,) 

then  h(tp)  [1  -  F(t)]  dt  {C«  -  Cp)  -  F(tp)  (C*  -  Cp)  - 

Cp  +  h(tp)  [*,Cp  -  »pC,l 

and  this  simplifies  to 

Now 

Maan  life  during  a  oyola 

Availability  A  (tp)  ■  - 

liq^eted  length  of  oyola 
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_ «,  tl  -  tflf)  <* _ 

It,  [1  -  #<y]  ■»  [«yi  t »,  [1  -  Fttpi  ♦  iA«  * 

_ L’’  p  -  ^yi  * _ 

U  -  «V1  -  K,  [«»p]  ♦  [1  -  «V]  * 

Setting  dJt  (tp) 

-  ■  0 

db 

we  get  (  U  -  F(mdt  +  R<IF<V  J  +  *^[1  -  F(V  ]  )  (1  -  r  (t^)  } 

■  «  [1  -  f(l)]  A  )  <  II  -  If  {t,)l  +  li,If(t,)I  -  M£U,)1  > 

Bimpli^ring  the  equation  above  we  get, 

QH~m\dt  +  *,[f(Vl  +  »,  -  l,If(t,>] 

»(V  [1  ‘  F(ff\dt  -  “  »»I»(t,)l  +  M'U,)1‘  * 

-  r(v  +i  > 

i«i£(s)i  -  <  } 

Cancelling  and  raairanging  the  terns  we  get 
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<£(v  I  [1  ]  c».  -  M }  -  {  »(tp>  [»f  -  RfCKt,)  j 

+  »,[»(*,)]  >  -  *,11  -  *<t,)] 
and  this  can  be  simplified  to 

(£(t,)  [  [1  -  flCOl  A  It*,  -  *,1)  -  (  '(tp)  [»,  -*,1  [1  -trtVl  - 

*,[1  -  F(t,)] 

Then 

—  f  ^  [1  -  F(i)  dt  -  F(tJ  -  — — 

1  -  F(,t^  Jo  ® 

[I  -  m  *  -  fti;>  - 
Rearranging  the  terms  we  have 
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APPENDIX  C.  DESCRIPTION  OF  IM8L  SUBROUTINE  DQDAG 


QDAQ/DQDAG  (Single/Double  precision) 

Purpose:  integrate  a  function  using  a  globally  adaptive  scheme  based  on 
Gauss'Kronrod  rules. 

Usage:  CALL  QDAG  (F,  A.  B.  ERRABS,  ERRREL,  IRULE.  RESULT, 
ERREST) 

Arguments 

F  •  User-supplied  FUNCTION  to  bo  integrated.  The  form  is  F(X),  where 
X  -  Independent  variable.  (Input) 

F  -  The  function  value.  (Output) 

F  must  be  declared  EXTERNAL  in  the  calllnp  program. 

A  •  Lower  limit  of  integration.  (Input) 

B  -  Upper  limit  of  Integration.  (Input) 

ERRABS  •  Absolute  accuracy  desired.  (Input) 

ERRREL  •  Relative  accuracy  desired.  (Input) 

IRULE  -  Choice  of  quadrature  njle.  (Input) 

A  Gauss-Kronrod  rule  is  used  with 
7-15  points  If  IRULE  -  1 
10-21  points  If  IRULE  -  2 
15  -  31  points  if  IRULE  -  3 
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20  -  41  points  if  IRULE  >  4 
25  •  51  points  if  IRULE  -  5 
30  <  61  points  if  IRULE  -  6 
IRULE  «  2  is  recommended  for  most  functions. 

If  the  function  has  a  peak  singularity  use  IRULE  >  1 
If  the  function  is  oscillatory  use  IRULE  «  6 
RESULT  >  Estimate  of  the  integral  from  A  to  B  of  F.  (Output) 

ERREST  •  Estimate  of  the  absolute  value  of  the  error.  (Output) 

Notes 

QDAQ  is  a  general-purpose  integrator  that  uses  a  globally  adaptive  scheme 
in  order  to  reduce  the  absolute  error.  It  subdivides  the  interval  [A,B]  and  uses  a 
(2k  1  )-point  Gauss-Kronrod  rule  to  estimate  the  integral  over  each  subinterval. 
The  error  for  each  subintarval  is  estimated  by  comparison  with  the  k-point  Gauss 
quadrature  rule.  The  subinterval  with  the  largest  estimated  error  is  then  bisected 
and  the  same  procedure  is  applied  to  both  halves.  The  bisection  process  is 
continued  until  either  the  error  criterion  Is  satisfied,  roundoff  error  Is  detected,  the 
subintervals  become  too  small,  or  the  maximum  number  of  subintervals  allowed 
is  reached. 
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APPENDIX  D.  DESCRIPTION  OF  IMSL  SUBROUTINE  DCSINT 

CSINT/DCSINT  (Single/Double  precision) 

Purpose:  Compute  the  cubic  spiine  interpoiant . 

Usage;  CALL  DCSiNT  (NDATA,  XDATA,  YDATA,  BREAK,  CSCOEF) 

Arguments 

NDATA  -  Number  of  data  points,  (input) 

NDATA  must  be  at  ieast  2. 

XDATA  •  Array  of  iength  NDATA  containing  the  data  point  abscissas. 
(Input) 

YDATA  -  Array  of  length  NDATA  containing  the  data  point  ordinates. 
(Input) 

BREAK  •  Array  of  length  NDATA  containing  the  breakpoints  for  the 
piecewise  cubic  representation.  (Output) 

CSCOEF  -  Matrix  of  size  4  by  NDATA  containing  the  local  coefficients  of 
the  cubic  pieces.  (Output) 

Notes 

DCSINT  computes  the  second  derivative  cubic  spline  interpoiant  to  a  set  of 
data  points  (x,,y,)  for  I  -  1,2 . NDATA  -  N.  The  breakpoints  of  the  spline  are  the 
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abscissas.  Endpoint  conditions  are  automaticaliy  determined  by  the  program. 
These  conditions  correspond  to  the  "not-a-knot”  condition,  which  requires  that  the 
third  derivative  of  the  spline  be  continuous  at  the  second  and  next-to-last 
breakpoint. 
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APPENDIX  E.  DESCRIPTION  OF  IMSL  FUNCTION  DCSITG 


C8ITQ/DC8ITQ  (Single/Double  precision) 

Purpose:  Evaiuate  the  integrai  of  a  cubic  spiine 

Usage:  CSITQ(A,  B,  NiNTV.  BREAK.  CSCOEF) 


Argument 

A  -  Lower  iimit  of  integration.  (Input) 

B  •  Upper  iimit  of  integration.  (Output) 

NINTV  •  Number  of  poiynomiai  pieces,  (input) 

BREAK  -  Array  of  iength  NiNTV-i-1  containing  the  breakpoints  for  the 

piecewise  cubic  representation,  (input) 

CSCOEF  •  Matrix  of  size  4  by  NINTV-t-l  containing  the  iocal 
coefficients  of  the  cubic  pieces,  (input) 


DCSiTQ  •  Vaiue  of  the  integral  of  the  spiine  from  A  to  B.  (Output) 


Notts 

DCSiTQ  evaiuates  the  integrai  of  a  cubic  spiine  over  an  intervai.  A  cubic 
spiine  is  a  piecewise  poiynomiai  of  order  4. 
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